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SYNOPSIS 


BAND STRUCTURE CALCULATION OP ELECTRON MOMENTUM 
DISTRIBUTIONS IN VANADIUM AND VANADIUM DEUTERI DE 

By 

ERISHNA GANDHI , X. R. 

Department of Physics 
Indian Institute of Technology 

Kanpur 208 016, India 

May 1278 

Recently the experimental determination of the electron 
momentum distributions (EMD) by the measurement of photon Compton 
scattering has become an important tool to study the electronic 
structure of solids. In order to make a detailed comparison 
between theory and experiment it is necessary to perform a band- 
structure calculation of the EMD. This is particularly important 
for the transition metals because of their interesting band- 
structure. The bandstructure and Permi surface of vanadium 
(which has partially filled 3d bands) are known to be sensitive 
to the choice of crystal potential, outer electron configuration 
ana exchange parameter. In the present work we have calculated 
the bandstructure and EMD (for photon Compton scattering) in 
vanadium for three different crystal potentials using the Hubbard- 
Mijnarends method. The results for the EMD are analyzed in terms 
of the bandstructure and Permi surface topology. Directional 

t 

Compton profiles (and their anisotropies) are calculated from the 
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EMI) obtained. The results are extended to calculate the 
directional Compton profiles for the vanadium deuteride in 
the framework of rigid band protonic model. 

In Chapter 1 a general introduction to the subject of 
EMjj is provided. The four important experimental methods to 
determine the EME, that is, photon Compton scattering, posi- 
tron annihilation, electron Compton scattering and (e, 2e) 
scattering are briefly described and compared with each 
other. A brief introduction to the theoretical calculation 
of the EME in metals is given. At the end the motivation 
behind the present work is explained. 

A review of the bandstructure and Fermi surface studies 
for vanadium reported in the literature is given in Chapter 2. 
dependence of these properties on the electron configuration 
and exchange parameter used in the crystal potential of vanadium 
is discussed. 

Chapter 3 provides a brief summary of Hubbard's fast 
approximate scheme used by us to calculate the bandstructure 
of vanadium. The extension of this method by Mijnarends to 
calculate the EMI) is also outlined. Necessary formulas and 
numerical parameters used by us are also presented. 

The results obtained by us for the bandstructure, Fermi 
surface and EME of vanadium for the three crystal potentials 
are presented in Chapter 4. The results for the EME are 
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discussed in terms of the bands true ture s Fermi surface and 
the group theoretical selection rule, fhe results for the 
EMI) along twelve p directions as well as the constant EMi) 
contours in the (100) and (110) planes are given. Possible 
ways of comparing theoretical and experimental EMjj are 
suggested. 


In Chapter 5 the results of the calculation of the 
directional Compton profiles ( p^) and their anisotropies 

arc described. I'hese results are also analyzed in terms of 
the bands tructure and Permi surface of vanadium. 


In Chapter 6 the results for the EMi) in pure vanadium 
are extended to obtain the EM n and the directional Compton 
profiles in vanadium deuteri des . 1'hese calculations were 
made in terms of a protonic model in which the dissolved 
deuterium is assumed to donate its electrons to the host 
metal. I'hese results are explained in terms of the band- 
structure of vanadium. A possible test to examine the 
electronic structure (and validity of the protonic model) of 
the transition metal hydrides is suggested. 



Chapter 1 


INTRODUCTION 


1.1 Introduction 

In recent years there has been considerable interest 
in the theoretical calculation and experimental determination 
of the momentum distributions of electrons in atoms, molecules 
and solids. At present there are four important experimental 
techniques for the determination of the electron moment-urn 
densities (EMD) and these are photon Compton scattering, 
positron annihilation, electron Compton scattering and (e,2e) 
electron scattering. These experimental techniques provide 
an indirect or direct information about the EMD, p(p), which 
is defined as 

p (?) = |x(?)l 2 (l.D 

where x(p) is the electron wavef unction in the momentum 
space. Usually the momentum wavef unction is obtained from 
the position space wavefunction i{)(r) through the well-known 
Dirac transformation 


x(p) = J'l'(r) exp (~ip*?) dr 


( 1 . 2 ) 
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The determination of the EMI), therefore, provides a test 
of the electron wave functions in atoms, molecules, metals and 
other solids and hence is of great interest. In the next 
section we shall discuss the technique of photon Compton 
scattering since the results of the present work can he 
directly verified by this technique. In §1.3 a brief outline 
of the other experimental methods to determine EI v Ii) will be 
given. The theoretical methods for the ca-lculation of the 
SMI) would be briefly reviewed in §1.4. Methods to compare 
experiment and theory are discussed in §1.5 while the moti- 
vation behind the present calculation would be discussed in 
§ 1 . 6 . 


1*2 Photon Compton Scattering 

Although it was known as early as in the 1930's that 
photon Compton scattering can provide information about the 
Mu, p(p), this technique lay neglected until 1965 when new 
experimental as well as theoretical developments led to a 
revival of research interest in this field. The interesting 
history of the development of photon Compton scattering 
technique is given by Cooper [1,2]. 

Photon Compton scattering is a name given to the inelas- 
tic scattering of a photon by an electron. In the case of a 
photon (wavelength X^) incident on a stationary, unbound 
electron the wavelength shift (ax) after Compton scattering 
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through an angle 1.1) is given by, 

4A = || sin 2 <f/2 (1.3) 

In a real solid the scattering electron is neither 
stationary (it has actually a finite linear momentum governed 
by the EMjJ) nor is it unbound. As long as the energy trans- 
ferred to the recoil electron is large compared to the binding 
energy of the electron, the effect of the finite nature of p 
(the electron momentum) is to modify Eq. (1.5) [1] as 

Ph r> 2 X -| 

AX = ££ sin */2 - (— sin */2)p z (1.4) 

where p^. is the projection of p on the scattering vector 

(taken to be along the z-*axis) . Examination of Eq. (1.4) 

shows that the intensity distribution l(X) of the Compton 

scattered line observed at an angle 4> will show a wavelength 

broadening instead of the monochromatic line as suggested by 

Eq. (1.5). The broadened lineshape l(x) v ersus x (Fig. 1.2) 

is usually termed as the Compton profile and the amount of 

broadening of a Compton profile is determined by the p -component 

z 

through the EMU [5,4]. 

The usefulness of the experimental Compton profiles for 
the determination of the EMh rests basically on the validity 
of the so-called 'impulse approximation (IA) 1 [3]. The 
central assumption in the IA is that the electron-photon 






6 


interaction time (involved in the Compton scattering) is short 
and the electron receives enough energjr to get excited into 
the continuum. This enables us to use a plane-wave descrip- 
tion for the final state of the electron. It has been shown 
by Platzman and Tzoar [5] that the IA is valid for free 
electrons while Eisenberger and platzman [6] demonstrated 
that it is also valid for the bound electrons provided that 
the electron’s recoil energy (%) is much greater than its 
binding energy (|E-g|). 

The cross section for Compton scattering (in the frame- 
work of the IA and using ft = c = 1) is given by 


dir 

dwafi 


ilSL -f p<5> 


a)2 


Th “l 


2 
2m 


-> 

LU 

m 


(1-5) 


ida 


where h~o| is' the elastic or Thomson cross section while 
a Th 


k^, cc-l (kg? “ 2 ^ l en °l e lb© momenta and energy of the incident 

°2* 


(scattered) photons (Pig. 1.5), and k = k-, - k OJ w = cu_, - cv 


^ 2 ’ w — x 

The cross section is thus related to the Compton profile 

J£(.q) 


(l) = J p(p) 8(q - ^~”-) dp 


1*1 


( 1 . 6 ) 


through the relation 


do 

dcodn 



f2 

“1 



Jg (q) 


(1.7) 
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where 


q = 


1*1 


\ 1*1 


( 1 . 8 ) 


and the quantity q is thus a measure of the electron momentum 
in the direction of momentum transfer. If the z-axis is 
chosen along i, the scattering vector, the expression for the 
Compton profile (CP) can he written as 


+ 00 


J£U> = // [p(p)] p =a d p x “p. 


(1.9) 


It is observed from Eq. (1.5) that the frequency of the 
scattered photon is shifted by the momentum transfer term 
(k 2 /2m) as well as by the Doppler shift term (k-p/m). The 
CP is thus obtained as a double integral [Eq. (1.9)] of the 
EMD, p (p) , in the momentum space. The relationship between 
the CP and the EMD assumes a simple form [Eq. (1.9)] in the 
framework of the IA and hence it is important to know the 
conditions under which the IA is valid. Investigation of this 
problem [6,7] has shown that the corrections to the IA are of 
the order of (E b /E r ) 2 where is the binding energy and E R 
the recoil energy of the electron. These considerations show 
that for vanadium (the metal studied in the present work) the 
IA is expected to be valid if gamma rays of energy E = 60 keV 
(e.g. from *"^Am) or higher are used. We have, therefore, 
performed the present calculations (Chapters 5 and 6) within 
the framework of IA. We shall now give a brief description 
of the experimental method. 
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The basic idea in the experimental measurement of a CP 

is to measure the energy spectrum of the photons Compton 

scattered at a. fixed angle § ^170° (Pig. l.l). Until recently 

X-rays (e.g. the 17.4 keV radiation o± the Mo— X— rays or 

the 22.2 keV radiation from the Ag-X a X-rays) were used as the 

primary incident radiation. The details of the X-ray Compton 

spectrometer are described by Cooper [lj. The use of X-rays, 

however, involved the following serious disadvantages: 

(i) The X-ray sources emitted radiation which has a doublet 

(i.e. K ) structure riding on a continuous bremss tra lilting 

cc j ~-a XI 

radiation. The primary incident radiation is thus not mono- 
chromatic. (ii) The wavelength analysis of the scattered 
radiation is usually done bjr employing a e B s2 9 B Bragg diffrac- 
tion measurement and this involves a slow point-by-point scan. 
This amounted to long counting periods typically of the order 
of several weeks. (iii) At the X-ray energies (typically < 

20 keV) the photoelectric absorption cross section is higher 
than the Compton scattering cross section and the useful 
(signal) counting rate is reduced. Secondly the validity of 
the IA which goes as (E^/E^) 2 becomes poor with X-rays if a 
scatterer with high Z is used. This restriction limited the 
use of X-rays to those elements which have Z < 20. 

These difficulties were considerably minimised when 
nuclear gamma lays were introduced as the primary incident 
radiation to study Compton profiles. Pelsteiner et al. [8] 
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2^1 

used the 59-57 keV gamma rays from the radioisotope r Am 
while Eisenberger and. lieea [9] used the 159-0 keV gamma rays 
from The Compton scattering experiments carried out 

nowadays use nuclear gamma rays from ^^Am or ^"^^Te or some 
other suitable radioisotope [5] and employ a solid state 
detector (e.g. Ce(li) detector) with the associated nuclear 
electronics to record the Compton profile data (Pig, 1.4). 

The main advantages of using nuclear gamma rays are their 
monochromatic nature, their high energy (which allows measure- 
ments with high Z scatterers within the IA) and higher counting 
rates (mainly because of the use of a solid state detector 
which allows a simultaneous recording of different scattered 
energies). The details of the gamma ray Compton scattering 
spectrometer, design of the experiment, processing of data and 
calculation of errors has been discussed in the literature 
[2,3,9]. 

It may be pointed out that recently the synchrotron 
radiation has been useu for the Compton profile studies [2,3]. 
However, these preliminary experiments only indicate the 
feasibility of such studies and more work is necessary before 
the synchrotron radiation becomes a routine source of the 
primary radiation. 

One important correction to be applied to the Compton 
profile data arises from the possibility of multiple scattering 


in the sample. The effect of multiple scattering on the Compton 
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profile is to reduce the intensity J(0) at q=G and to broaden 
the profile [3j. A correction for the multiple scattering 
using a Monte Carlo method has been suggested [3,10], Accurate 
Compton profile measurements require a choice of optimum 
thickness of the sample and proper correction to the data 
collected, 

1.3 Som e O ther Experimental Met h o os fo r the Determination 

o f Md 

1.3.1 I ntro duction 

In this section we shall briefly outline three other 
important experimental methods employed for the determination 
of the SMI. These methods have been described and compared in 
detail in the literature [11-15]. Our objective here is to 
provide a- brief introduction to these methods for the sake of 
completeness . 

1.3.2 Po sitr on Annih ilation 

Upon entering a material medium positrons lose their 
energy by ionization ana excitation and become thermalised in 
a time ('vlO sec) which is 3. ess than the lifetime of posi- 
trons sec or longer) . let us consider the case of 

metals in which such thermalised positrons annihilate with 
the electrons mostly in a singlet (^S ) state with a lifetime 
% ^ 10“ 10 sec or higher, buch an annihilation is accompanied 
by the emission of two photons. The lifetime of the positron 




rv 
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together with the direction of emission (or momentum) and 
energies of the annihilation photons carry information about 
the positron-electron state prior to annihilation. The 
finite momentum (in the centre -of -mass frame) of the positron- 
electron pair arises almost entirely out of the electron 
momentum (the positron being thermalised gives k + * 0) and 
gives rise to the photon pair momentum jT. The finite value 
of p manifests itself in two ways. The transverse component 
Pj (perpendicular to the direction of emission of the photons) 
gives rise to a deviation from collinearity in the path of the 
annihilation photons. If p^ = 0, the two photons are emitted 
exactly at 180° apart, but if p x / 0 the two photons deviate 
from collinearity such that 0 = p /me and 0 = p /me 

^ ^ V •J 

(fig. 1,5). for the values of p encountered, these angles 

and 0^ are of the order of a milliraaian (one milliradian 

-3 

corresponds to a momentum of me x 10 or 0.137 atomic momentum 
units). 


The longitudinal component p j j (parallel to the direction 
of emission of the photons) is responsible for a change in the 
energy of the annihilation radiation according to the relation 


Y 


= me 2 + 

8m 




fit 


( 1 . 10 ) 


Two experimental methods based on these two effects are used 
to study EMD by the positron annihilation method. In the 
first method, called the angular correlation of positron 
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annihilation radiation (ACPAB.) , one determines the angular 

distribution of the coincident counting rate IT as a function 

of the angles 0 and 8 . If e__ is left undetermined (i.e., 

z y y 

the detectors have a long dimension along the v-axis in the 

form of a . long slit so that e_ T » 0 ) one measures the H(6 ) 

y z z 

curve alone and such a geometry is called the long slit 
geometry [12 j. In an actual experiment the source of positrons 
is several millicuries of the radioisotope Na or Cu (which 
emits positrons of low energy <400 keY) . fhe measured 
coincidence counting rate ll(p = me e ) is related to the 
EMh, p(p) , through the relation 

+ oo 

*(P Z ) = // P(p) dp x dp (1.11) 

— CO v 

which is same as the relation [Eq. (1.9)] between the CP, J(q) 
and the EI-ID, p(p). The momentum resolution obtained by the 
long-slit geometry is of the order of 0.07 a.u. whereas the 
corresponding figure for photon Compton scattering is 0.4 a.u. 

Recently the Brandeis group [3,16] has developed a multi- 
counter two-dimensional angular correlation apparatus with 
which lT(p , p ) curves can be measured. 'Their apparatus 

y z 

consists of a (32 x 32) matrix of rectangular llal(Tl) counters. 
Such a two-dimensional curve l\T(p , p ) measures a line integral 

y z 

of the BMP through the relation 

= / p(p) d P x 


lT(Py,P z ) 


( 1 . 12 ) 



15 


As compared to the Ion/;, slit geometries the two-dimensional 
AC PAR curves provide a greater wealth of inf ormation. Various 
other geometries have "been used in designing an ACPAR measure- 
ment ana these are describee in the literature [11,12]. 

Although the positron annihilation experiments ofler a 
better momentum resolution when compare a to photon Compton 
scattering and a choice of different geometries, they suffer 
from two limitations, firstly, the crystalline samples used 
for positron annihilation measurements must be free from 
defects otherwise these defects can trap positrons and alter 
the angular correlation curves [11,17]. Secondly, the theore- 
tical interpretation of the positron results requires a 
knowledge of the positron wavefunction and positron-many 
electron correlation effects. In spite of these limitations 
the experimental data obtained by the positron annihilation 
studies of transition metals have compared well with the 
results of bandstructure calculation of the EMI) [18]. 

1.5.5 Electron Compton Scattering 

Inelastic scattering of fast (20-50 keV) electrons is 
termed as electron Compton scattering or electron impact 
spectroscopy. Although such measurements were carried out 
in the 1950 's this technique lay dormant until about 10 years 
ago when the new developments in the theory and experiment 
helped to revive this field [15,19]. In principle. 
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the experimental method, to study the electron Compton scat- 
tering from gaseous samples consists of shooting a collimated 
electron beam (incident energies of the order of 20-50 keV) 
at the gaseous sample which is injected in the form of a 
narrow jot in an evacuated chamber (fig. 1.6). The scattered 
electrons are energy-analyzed by a Mollenstedt velocity 
analyzer and counted by a Si solid state detector to obtain 
the l(E) versus S curve (or the electron Compton profile). 

'Hie scattering angle can bo changed from 0° to 140 ° by rotation. 
In such a measurement of the BKu, a resolution of 2-5 eV 
(equivalent to a momentum resolution of 0.01-0.03 a.u.) is 
usually obtained. In addition to such good momentum resolu- 
tion, the counting rates are quite high and the final accuracy 
is of the order of 1 % or so. 

There are some theoretical assumptions and corrections 
which restrict the theoretical interpretation of the data 
obtained by electron Compton scattering [15,19]. The first 
basic approximation consists in applying the Born approxima- 
tion to describe the inelastic scattering process. This 
approximation is certainly valid for small atoms and molecules. 
Secondly, the binary encounter approximation is used to 
calculate the differential cross section. The latter approxi- 
mation assumes the final state of the electron to be taken as 
a plane wave and is analogous to the IA used in the photon 
Compton scattering. Next the cross section calcuilated using 
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the above approximations has to be corrected, for (a) the 
relativistic effects and (b) the exchange and interference 
between the incident and scattered electrons. 


It has been shown that within the above framework the 
Lifferential cross section can be written [13] as 


2k s [1 


(1 - |3 2 )V2c 2 ] 


afidli 


k ± (l— f3 2 ) [K 2 - S 2 /4 c 2 ] 


E ex J(<l) (1 ‘ 13) 


where k- and k_ describe the wavevectors of the incident 
1 s 

and scattered electrons, lii is the energy loss on scattering* 

■f 

E is the momentum transfer, K = k. - k and q = [E(l-E/4c) -K 2 ] , 


ihe function E takes into account the exchange and inter- 

G jL 

ference corrections and 



while J(q) is the Compton profile which, for the case of 
isotropic EMD < P(p) > , assumes the following form 


J(q) = 2 % f <p(p)> p dp (1.15) 

fq/ 

Uptill now the usefulness of the electron Com.pt on scat- 
tei*ing technique for the study of EMD has been restricted to 
the g-aseous samples of the elements like H and He. In terms 
of the momentum resolution, statistical accuracy and counting 
rates available, this method is definitely superior to other 
methods for studying gaseous samples of low Z. However, this 
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method has serious experimental ana theoretical difficulties 
for extending it to heavier gases (e.g. Kr, Xe etc.) or to 
the bulk (or solid) samples (where the multiple scattering 
of electrons can be a serious disadvantage) . 

1.5 . 4 (e,2e) Method 

In the photon (and electron) Compton scattering experi- 
ments described earlier the information carried by the recoil 
electron is lost because only the energy of the scattered 
photon (or electron) is analyzed. In the so-called (e,2e) 
method the energy and momentum of the incident electron as 
well as the knocked out electron (both detected in coincidence) 
are determined. The experimental details and the theoretical 
framework of the (e,2e) method are reviewed by McCarthy and 
Weigold [15]. The basic geometry of the (e,2e) method is 
shown in Tig. 1.7. The incident electron having momentum k Q 
and energy E Q knocks out an electron from the target. The 
two electrons (scattered and knocked out) coming out with 
momenta k^ and kg and having energies E^ and Eg respectively 
are detected in coincidence and the coincidence counting rate 
can be measured as a function of the parameters E Q , E^, E^, 6 
and <t> (Pig. 1.7) • The results of the measurements are 
analyzed in the framework of the 'distorted wave off -shell 
impulse approximation' [15] in which one uses slightly 
distorted piano waves to describe the incoming and the two 
outgoing electrons and the distortion is calculated within 
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the eikonal approximation. The differential cross section 
is then ca.lculated [13] from 

= vtt-, 2 / 2me 2 \ 2nn r 

d dHg dE^ nl ^ ^2 ' [exp( 2%v,) -lj 

{ j k-rk ; | 4 | k-k' | 2 |£+k'| 2 

x cos [p In (i|t|Lif)]> | X (E)| 2 (1.16) 

|k-k ! | ^ n± 

where b ^ is the spectroscopic factor and 

N = ex P HVo + Va + Vb )R 3^ 

r] = me 2 /h 2 k' 

k = (E q -E), k'= (K a -E b ) 

K = K 0 - K a - K b 

t 2 ± = (1 + p + iY)k ± , 1 = 0, k, B 

and the subscripts 0, k, B stand for the incoming and two 
outgoing electrons. 

The cross section (1.16) is thus seen to be directly 
proportional to the EMu itself and not to a single (or double) 
integral of the EMU as in the previous methods. Apart from 
this major advantage, another advantage of this technique is 
that one can determine the EMx> , P n ]_(p)? cue to each individual 
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orbital (nl) . In spite of these advantages, this method has 
so far found limited applications because of the following 
experimental difficulties; (a) only gaseous samples can be 
stuaied ana therefore only isotropic (or radial) momentum 
distributions can be determined. (b) Owing to the low counting 
rates from the core orbitals they cannot be studiea and only 
the outer electron states can be examined. In view of these 
difficulties the (e,2e) method has been so far applied only to 
gases and thin films. 

1.2.5 Comparison of the . Vari o us Experimental Methods for the 
.Determ i nation of the~ BMP 

The main advantages as well as the limitations of each 
of the above four methods have already been pointed out in the 
previous discussion. In Table 1.1 we provide a general 
comparison of the four methoas described above. It is now 
well recognised that in spite of its moderate resolution, the 
technique of photon Compton scattering is the best single 
method suited to study the effect of electron wavefunctions 
on the EMI in a variety of substances. ’The positron annihila- 
tion method offers the advantages of superior momentum 
resolution, choice of geometry, possibility of measuring the 
two-dimensional ACPAE curves etc. but it is restricted in the 
theoretical interpretation available. The positron annihila- 
tion technique is more suited as a method to investigate the 
Fermi surface of metals and metal alloys. The other two 



Table 1,1: Comparison of different methods 
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methods, electron Compton scattering and (e,2e) scattering- 
have not yet proven their complete potential hut it is 
doubtful whether they can ever he used tc study the bulk 
(or solid) samples, ihe best approach seems to study a 
particular sample both by photon Compton scattering and 
positron annihilation and combine these results to determine 
the EMjJ in the system. 

1 Calculation of the electron Mo me n t urn D e ns ity 
Choice of the Electr o n lavef unction 

•Hie first step in the calculation of the Ebn, P(p), is to 
choose the e3.ectron wavefunction appropriate for the system 
under study. As a crude approximation one could use a free- 
electrcn model or free atom wavefunctions [3j. However, it 
has now been demonstrated [3] that for metals it is essential 
tc make a bands tructure calculation employing one of the 
standard (e.g. 01 W, APil, KKH or other) methods, for other 
solids (e.g. liP) one may use a methoc rJ<- sed on the linear 
combination of atomic orbitals [ICAO]. In the case of th.e EMI) 
fox positron annihilation iu is requii-ed to know the positron 
wavefunction as well. The positron wave function is usually 
determined by the ¥igner -Seitz method or the plane-wave 
expansion method or by the positron pseudopotential method. 
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1. 4.2 Electron-electron Correlation 

Usually the bandstrueture calculations of the electron 
wavef unctions are performed within the independent particle 
model. In this model the electron-electron interactions are 
reuuced to single particle potentials, ihe electron-electron 
interaction arising from the antisymmetry requirement of the 
electron wavefunction is termed the exchange correlation part 
and that arising from the Coulomb ic repulsion between the 
electrons is called the pure electron-electron correlation. 

It has been shown [5] that the electron correlations can 
significantly affect the EMl) and the CP curves for the simple 
[e.g. Li, Na, Al, ...] metals. In the case of transition 
metals one usually neglects the pure electron -electron cor- 
relation effect. In the present calculation we have taken 
into account the exchange correlation explicitly and neglected 
the pure electron correlation. 

1.4.5 Theoretical Models 

In the following sections we shall briefly outline the 
theoi'etical methods for calculating the EKE for a free- 
electron gas , free atoms and metallic solids. 

( A ) Free electron gas 

The EMI), p(p) , for a non-interacting electron gas is 
isotropic in nature (i.e. p(p) = <p(p)>) ana at I=0 o K it 

assumes a constant value up to p <_ p ^ (where Pp is the Fermi 
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momentum) and it drops to zero in a discontinue -us manner 
(Pig. 1.8a). The corresponding Compton profile is parabolic 
in nature and is given by 

J(P Z ) = 2m (p* - p|) for p z < p p 

= 0 for p z > pp (1.17) 

The effect of taking the electron-electron correlations 
into account (i.e. interacting electron gas) is to transfer 
a part of the occupied states (lying just below p=p^) to the 
higher momentum region (just above p = p™) as shown in 
Pig. 1.8b. This results in a smoothening of the disconti- 
nuity in the EMn, p(p) , at p = p^, and an appearance of a 
tail beyond it. 

( B) The free atom 

To gain an insight into the EMu contributions by atomic 
orbitals we have compared in Pig. 1.9 the radial charge dis- 
tribution in argon in the position (r) space and the EMI) in 
the momentum space. The reciprocality between the position 
and momentum space (i.e. a narrow function in the r-space 
becomes a broad function in the p-space ana vice versa) is 
brought out by these distributions. Thus, for example, the 
EMI for the M-shell electrons appears sharp although its 
position charge density appears diffuse. 
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While comparing the experimental Compton profile data 
with the band theoretical results one needs to know the 
contribution to the El'in from the inner core orbitals which are 
free -atom-like, ihis enables one to compare only the 'band- 
electron' profiles. The El In due to the atomic orbital (nl) 
is assumed ro be isotropic ana is usually calculated from the 
f ormula 

P n i(p) = (§) f / dr r 2 P lll (r)5 1 (pr)| 2 (1.18) 

where (P^/ 1 ') are the atomic wavef unc t ions (say, for example, 
those calculated by Herman and Slcillman [20]). Recently, 

Biggs et al. [21] have published Tables of theoretical 
(Hartree-Eock) Compton profiles for the atomic orbitals of 
various elements. The various aspects of the atomic Compton 
profiles are discussed by Ifenaelsohn and Sfaith [22]. 

( c) Solids 

In a solid, the Compton profile is sensitive to the EMh 
(and the wavefunctions) of the ou.ter electrons and thus 
provides a valuable technique to study the solids where these 
outer (or conduction) electrons are responsible for the 
bonding. The discussion of the EMU in solids can be divided 
according to the type of solid (e.g. metal, alloy, metal 
compound, covalent crystal, ionic crystal etc.). In an 
excellent review article Ber^gren et al.j_23] have reviewed 
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the theoretical and experimental studies of the Compton 
profiles of solids. Their review neatly summarises the situa- 
tion regarding the 3d transition metals, the area of our 
present interest, file free-electron gas or a free-atom model 
is known to give a poor description of the EMi and CP of 3d 
transition metals primarily because of the complicated effects 
of the (4-s+3d) electrons. The disagreement occurs mainly in 
the low momentum region where the major contributions come from 
the outermost electrons. 

Berggren [24] in 1972 applied the renormalized-f ree-atom 
(REA) model [25] for the calculation of the EMI) and CP in 7. 

In the I: PA model one proceeds as follows. The free-atom 
(e.g. Hartree-Bock) wavef unctions are truncated at the radius 
r TT g of the II igner -Seitz (¥S) sphere and then renormalized to 
unity within the PS sphere. Such a procedure takes into 
account the fact that in a solid the electron wavef unctions 
are contained in a cell. It has been shown [25] that such 
RPA wavef unctions provide a good estimate of some of the band- 
structure properties of the 3d transition metals. The exten- 
sion of the PlPA model for the calculation of the EMD and CP 
is described by Berggren [24] and it has been shown to give 
satisfactory agreement with experiment for the transition 
metals like 7 [23] s Ti [23], Pe [23], Ei [26] ana Cu [26], 

In spite of its success the main limitation of the RFA model 
is that it can provide onljr a gross description of the CP for 
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the polycrystalline samples hut cannot provide a proper des- 
cription of the finer features like the anisotropy in the 
directional profiles (from single crystals) „ In the transition 
metals ( parti cularly for a metal like 7 for which the 3d band 
is incomple tely filled) the topology of the fermi surface and 
the complicated nature of the 3d vavefunctions is reflected 
in the anisotropy in the UHL and OP. It is, therefore, 
essential to carry out the calculations of the Mu with the 
help of some standard bands tructure methods. 

As mentioned above, the bands tructure calculations of 
the EHjJ in a solid are usually carried out in the framework 
of independent particle model (I®) in which the expression 
for the Blip p(p) is expressed as 


occ . 

p(p) = I. f(k,j) Ia - (k, p) 1 2 o( p-k-K - ) (1.19) 

j ,1c, i 3 1 

In Eq. (1.19) f(k,j) is the Ecrmi-Pirac distribution function 

* 4 * y 

for an electron in state k and band j such that at T = 0 K, 
f(k,j) = 1 if the state (k,j) is occupied ana f(k,j) = 0 if 
the state (k,.j) is unoccupied. The set {E^> is the set of 
reciprocal .lattice vectors and 




p) = 


/ 

cell 


-ip 


- r 


*k,f r) d? 


( 1 . 20 ) 


In Eq. (l. 19) the summations over j , k and take into 
account the sum over occupied states and the normal (K^=0 

4- -> Hh -4* 

and p = k) as well as the Umklapp (K^o, p = k+K^) processes 
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thus bringing out the difference between the real physical 
momentum p and the crystal momentum It. 

An electron in a state it, thus contributes not only 

-V 4 1 *+■ 

at p = k but also at other p = k+K^ with an amplitude 
A.(k, p). The importance of the Umpklapp (K^ 0) coefficients 
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A^(k,p) (also termea in the literature as high -momentum- 
components -HMC - or the high— Fourier -components HFC) can be 
illustrated by following Berlto and plaskett [27]. 

For the sake of illustration let us consider the simple 
case of nearly -free -electrons (l-JFB) in the one -dimensional 

-y . -K 

model. A comparison of the k- and p--space distributions P(k} 
and p(p) is shown in Fig. 1.10 for the following three cases 2 
(a) Fermi surface lying in the first zone, (b) 'So Fermi 
surface 5 p(k) distribution touching the boundary of the 
JBrillouin zone (BZ) and (c) Fermi surface in the second band 
(shown in the extended zone scheme). It is observed from 

. *■>* V 

Fig. 1.10 that the breaks in the SHU, P(p), occur not only 

-y ->* 

in the first zone (lh=0) but also in the higher (K^O) zones 
although with a reduced Umklapp amplitude A-(k,p). In the 

O 

case of a transition metal the EMI (ana hence the CP) is 
therefore determined by the Fermi surface topology and the 

-y 

k- (and p-) dependence of A.(k,p). There is also an important 

j 

effect due to the symmetry selection rule [28] which we shall 
discuss in Chapter 4. These considerations show that in 
order to understand the basis of the Compton profiles in a 
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metal, the theoretical results of the EILu have to be carefully 
examinee. in terms of its banes true ture and Fermi surface. 

1 . 5 Comp arison Between Theory and experiment 

There are two possible approaches to compare the experi- 
mentally measured Compton profiles with theoretical results. 

In the first method one starts with the theoretical EMjJ s p(p) 
and obtains the Compton profile J(p ) from the following 
relation; 

J(pj = JJ p(p) dp (1.21) 

ihe integration involved in (1.21) requires a knowledge of the 
EMh, P(p) , at several p-points in the momentum space and 
therefore involves a large computational effort. An approxi- 
mate but faster way of computing the integral (1.21) is to 
expand p(p) in terms of ICubie harmonics. Vie have followed the 
latter approach to calculate J(p^) and this method would be 
described in Chapter 5. 

Another method consists of unfolding or inverting the 
EK1), p(p)j from the experimentally measured set of the 
directional profiles (p ) (i.e. for a set of {hkll single 

crystal directions) through Eq. (1.21) . Two such methods 
have been suggested in the literature by Mijnarends [29] and 
Mueller [30]. The results of the 'unfolding' analysis lead 
to the p(p) distributions which can then be compared with the 
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theoretical distributions. ihese methods of the reconstruc- 
tion of the EKD ere reviewed by I-Iijnareiids [51]-. 

1 . 6 Motivation Behind the Present Work 

It was pointed out earlier that the technique of photon 
Compton scattering is well-suited for testing the electron 
wavex unctions via the EMu. For ferromagnetic iron liijnarends 
[32] has compared the EMl) obtained by unfolding the experimental 
positron annihilation data with the theoretical EMjj, calculated 
from the bands tructure of Wood [33] and Vakoh ana Yamashita [34]. 
His results demonstrated that the experimental data support 
the bands tructure of Pe by wood [33] . 

Similar tests can be carried out on other metals, espe- 
cialljr those transition metals for which the Fermi level lies 
among the a bands . The character of the d wavefunctions , 
which is mainly responsible for the anisotropy in the EMI 
(and CP) can provide a sensitive tost for distinguishing 
among different bandstructures . One such transition metal is 
vanadium in which the d bands are incompletely filled- The 
results of bandstructure and Fermi surface of V (described 
in detail in Chapter 2) show interesting changes with respect 
to the different crystal potentials used in the literature. 

‘The 5 band electrons from the (3d) il (4s)^“ ri configuration fill 
up the entire first, most of the second and about half of the 
chird conduction band in V. As a result the Fermi surface (FS) 
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of V shows a complicated structure and the IS topology (and the 
bands tructuro) is sensitive to the crystal potential and the 
bandstructure method used in the calculation. Vie, therefore, 
thought it interesting to undertake a banastructure calculation 
of the SMi) in V and to examine the effect on the EICD by chang- 
ing the crystal potentials Although the Compton profiles of 
metallic V have been theoretically calculated by bakoh and 
Yamashita [35]? their work does not give any details of the 
behaviour of the EMi) itself. In the present work we have 
tried to examine the EMD in terms of the bands tructure ana FS 
effects. In other words we treat the EMU as one more ground- 
state property (similar to electronic specific heat, suscep- 
tibility etc.) predicted from the bandstructure calculations . 

Ue have also extended the results for the EMD to obtain the 
directional Compton profiles and their anisotropies. 

The experimental Compton profile data on poly crystalline 

li was analyzed by Berggren et al. [23] in terms of the UFA 

model and the results were interpreted to show that the 

9 *2 

electron configuration of Ti is closer to (3d) “(4s)^ rather 
than (3d)^(4s)l. One of the early calculations by Mattheiss 
[36] has shown that the bandstructure of V changes considerably 
if the electron configuration is changed from (3d)^(4s)^ to 
(3d)^(4s)^. In order to examine the changes brought in the 
EMn (as well as in the directional CP curves) by the changes 
in the outer electron configuration, we have calculated the 
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EMjj using two crystal potentials one for the electron configu- 
ration (3d) 4s) 2 and other for (3d)^(4s) 1 . A third type of 
potential was also used. 'These details would he discussed in 
the following chapters. 

Recognising the interest in the field of EMJJ, the 
International Union of Crystallography .formed a Special 
Commission on Electron Charge, Spin and Momentum Density in 
1972. A special project was initiated by this Commission in 
1975 to determine the accuracy" of tlio theoretical Compton 
profiles and X-ray scattering factors for metallic V. The main 
object in this project has been to compare the results of 
various theoretical methods and to determine their relative 
accuracy thus providing the experiments. lists a basis for 
comparison.' 

Our research group joined this International project and 
the present work was •undertaken with this object in mind. The 
standard lattice potential of metallic Y to be used for this 
project was prepared by Mi jnarends [37] for the electron 
configuration (3&)^(4s)^. \!e decided to extend the scope of 

this project by calculating the EI-Lu and CP in Y for two other 
crystal potentials besides this standard potential, 

Recently the electron structure of the hydrides (and 
aeuterides) of the transition metals have invoked considera- 
ble interest. In the metal hydrides the electronic states of 
the outer or conduction electrons are modified by the dissolved 
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hydrogen (or deuterium) . The technique of photon Compton 
scattering provides a sensitive method to probe the conduction 
electron states ( via the EED) of a metal hydride and recently 
the Compton profiles of hydrides and deuteriues of niobium 
[i.e. KfbH and Nbi) ] have been investigated in detail [38]. 

We have therefore found it interesting to extend our theore- 
tical results on the EME in pure V to VE, r to calculate the EME 
and directional Compton profiles in VE„ . It is well known 

-A. 

that the group VB metals V and Kb have similar bandstructure 
and Eermi surface topologies. Hence a calculation of the EME 
in Vi) a.nd its comparison with HbE^, as well as with experi- 

-A. _X_ 

ment, can provide valuable information about the electronic 
structure of transition metal hydrides. 

Keeping these objectives in mind we have calculated the 
bandstructure and EMI) in V applying Hubbard's fast approximation 
scheme for three different types of the crystal potential. 
Chapter 2 provides a review of the bandstructure and Fermi 
surface of V as reported in the literature and summarises the 
motivation behind our work. A brief description of the 
Hubbard -Mi jna rends method of ca icula ting the bandstructure 
and EME is given in Chapter 3. The theoretical results of our 
calculations are presented in Chapter 4 where these results 
for the EMI) are analyzed in terms of the bands tructure and FS. 
Directional Compton profiles and their anisotropies are 
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obtained from the EMI and are presented in Chapter 5 • The 
extension of our results to vanadium aeuteride ( Vv) and the 

2L 

calculation of the directional Compton profiles in 1/1 are 

Ji. 

discussed in Chapter 6. 
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Chapter 2 


REVIEW OP EA.Fi/S IRUC TIRE AMD FEEI'iI SURFACE 
STUDIES OP VAFADITM 


2.1 Introduction 

There has been considerable interest in the electronic 
structure of the "vanadium group (i.e. the VB isoelectronic 
group consisting of vanadium , niobium and tantalum) of metals. 
In all these three metals the outer (or conduction) electrons 
fill up the outermost d-states (e.g. 3d for V) only partially 
and the Fermi level straddles the d bands thus resulting in 
a complicated topology of the Fermi surface (PS). This has 
resulted in a number of theoretical calculations of the 
bandstructure and experimental investigations of the Fermi 
surface of the group VB metals. Before we review these, 
theoretical and experimental studies for V (the metal of our 
interest) we shall outline some features of the bandstructure 
calculations which are of interest in the discussion of our 
results (Chapter 4). 

The calculation of the energy levels of the electrons 
in a solid is based on the basic assumption (involving the 
Born-Oppenheimer approximation) that the electronic and 
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nuclear notions can be describee separately and. that the 
effect of the nuclea.r notion on the electrons is sma.ll and 
negligible . As atoms are brought together to form a crystal- 
line solids the potential assumes the periodicity of the 
crystal lattice . The calculations are carried out in the 
independent particle model (IM) in which one limits the 
ca.lcula.tion to the states of a single electron in a rigid, 
infinite ana periodic lattice. The results for the many- 
electron system, are then derived by treating the electrons 
as a system of non-interacting, particles in this periodic 
potential and applying the i'emi-nirac statistics to 
accommodate them in these states. Many bulk properties of 
the solids can be fairly successfully described in the 
framework of the above model. 

A critical discussion of the various approximate 
methods used to solve the bandstructure problem is beyond 
the scope of this work. Excellent review article [1] and 
texts [2-5] exist in the literature on this subject. In the 
following we shall outline orQgr those features which are 
relevant for the discussion of our results. In this work, 
unless otherwise stated we shall follow the Hartree system 
of atomic units in which h = 1, m = 1, e=l and c = 137.036. 

file Hartree-Eock -Slater equations 
H #(r) - [-■£-+ V c (r) Y x ( r) ] *(r) 

= [- k + V(r)] ♦(?) = S *(r) 


( 2 . 1 ) 



44 


provide a basis for the solution of the bands true ture problem. 
In Eq. (2.1) V_ represents the exchange potential arising 
out of the fact that electrons have spin, are indistinguishable 
and obey the lauLi exclusion principle. On the other hand, 

V c (Eq. 2.1) is the non-exchange potential experienced by 
the electron in the Coulomb field of all ions and all other 
electrons . 

The next problem in the band theory is how to solve the 
one-electron equation. This problem can be divided into two 
parts s (a) the mathematical method to be used and (b) the 
one-electron potential, V(r) in Eq. (2.1), to be used. Any 
one of the methods (e.g. the matrix method, the variational 
formulation or Green's function method) can be taken as the 
mathematical method to be used. In the present work we have 
used Hubbard's fast approximation scheme (which is based on 
the EKRZ method) and it is described in the next chapter. 

After making the choice of the mathematical method to 
be used, one has to make the important decision of selecting 
the initial form of the potential V(r) involved in Eq. (2.1). 

By employing an iterative procedure one can modify V(r) in 
each iteration so that a satisfactory self-consistency (SC) 
is achieved between the V(r) and the calculated wavef unctions . 
Alternately one may perform the calculation non-self-consis- 
tently (NSC) . ¥e may expect the following behaviour of 
?(r) as atoms are incorporated into a solid, (a) The potential 
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V(r) remains spherically symmetric [i.e. Y(r) = Y(r)] in 
the region of the ionic cores, (b) The potential Y(r) remains 
fairly constant in the interstitial regions outside the 
ionic cores. These ideas are incorporated in the muffin- 
tin (MT) type of potential which assumes Y(r) = V(r) inside 
spheres of radius centered on each nucleus, and Y(r) = Y q 
outside these spheres. Such a design of the MT potential 
makes it possible to solve the one-electron equation exactly 
and accurately. Within the inscribed sphere the equation 
gives separable solutions in polar coordinates while in the 
interstitial region (where Y(r) = V ) plane-wave solutions 
are used. The constant potential V Q is usually taken as the 
average value of the potential between the boundary of the 
inscribed (MT) sphere and the V/igner-Seitz cell. The radius r^ 
of the inscribed sphere is chosen arbitrarily so that the 
spheres do not overlap. 

The construction of the muffin-tin crystal potential 
usually follows the method suggested by Mattheiss [6], The 
Coulomb contribution, V (r) (in Bq. 2.1), is obtained as a 
superimposition of neutral atom Coulomb potentials (arising 
from the neighbouring atoms) expanded about the origin using 
155w din's alpha function expansion method [7]. Atomic orbitals 
and atomic potentials obtained from Hartree-Fock-Slater self- 
consistent field equations (e.g. Herman and Skillman [8]) 
are usually used. The super-imposition is done in two stages $ 
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the Coulomb and exchange contributions are considered sepa- 
rately. The exchange potential term V x (r) is usually calcu- 
lated in the free-electron-exchange approximation of Slater 
[9,2,3] who assumed that the exchange term for the electrons 
in the potential field of crystal is similar to the one for 
electron gas and depends only on the local electronic charge 
density p(r) . Ihe averaged exchange potential is obtained 
from 

V z (r) = -6 a [(3/8*) |p(r)|] 1/3 (2.2) 

Rydbergs 

There is some arbitrariness about the proportionality 
constant cc 5 which can be taken as a= 1 (Slater's full exchange 
[9]) or a= 2/3 (Gaspar [10], Eohn and Sham [11]) or cc can be 
varied between these limits to improve the quality of the 
bands tructure results. 

A particular bands tructure calculation is, therefore, 
generally characterized by (i) the particular mathematical 
method used (e.g. KKR, AP¥ etc.), (ii) the self-consistency 
or the non-self-consistency requirements for the potential 
V(r), (iii) the model for the crystal potential (e.g. muffin- 
tin) , (iv) the electron configuration used for the outer 
(band) electrons in constructing the crystal potential, 

(v) the value of the exchange parameter a used in Eq. (2.2) 
and so on. 
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The aims of the final calculation are as follows : 

(a) To calculate the energy relations E(k) versus k along 
different k directions of the single crystals of the 
solid. 

(b) To determine the iermi energy Bp through the calculation 
of the dens it y-of -states curve n(E)dE and use this result 
to construct the constant energy surface , B(k) = Ep 
(i.e. Fermi surface) in the $ -space. 

(c) To determine the wave function .(r) describing the 

electron belonging to the jth band in a particular 
k-s tate. 

( d) To combine the results of (a), (b) and (c) above to 
determine various physical bulk properties of the solid. 

Nowadays the bandstructure problem is meant to compre- 
hensively cover the above calculations. It is only recently 
that these calculations are being extended to theoretically 
determine the electron momentum distributions, p(p) (Eq. 1.18)* 
as one more ground state property predicted from the band- 
structure results. 

The transition metals constitute a large group of metals 
in the periodic table and there is considerable interest in 
their electronic structure and other physical properties in 
view of their technological applications. In the iron series 
(or 3d) transition metals the outer 4s electrons fomn a broad 
s band which behaves like a nearly -free-electron (EFE) conduction 




Hybridisation or tlie a arid d bands 
for a fee transition metal. 
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band, in an ordinary metal. Fhe 3d electrons form a narrow 
3d-band because of the overlap with the neighbouring atoms. 

The narrow d band lies within the s-p (or conduction band) and 
hybridizes with it when they cross (Pig. 2.1) . The relative 
positions of these bands (and therefore the bands tructure, 
Fermi surface and the electron wave! unctions) in the transi- 
tion metals are sensitive to the choice of bandstructure 
method, crystal potential, outer electron configuration, 
exchange parameter a and so on. In the present work we have 
investigated the sensitivity of the EtfD, P(p) in vanadium to 
the choice of crystal potential (including the choice of 
electron configuration and the exchange parameter cc). 

In the next two sections we shall review the results 
obtained by other workers on the bandstructure and Fermi 
surface of vanadium. This background would be useful in dis- 
cussing the results of our calculations. 

2 . 2 purv ey o f the Bands tructure Calc ulation s for Vanadium 

The earliest bandstructure calculation on vanadium was 
made by hattheiss [12] who reported the sy sterna tics of the 
bands t rue tures (along a line of symmetry) for the iron-series 
transition metals using a non-self -consistent (NSC) augmented- 
pl ana «wave (AW) method. His results demonstrated that for V 
the bandstructure (in particular the location of the. d bands) 
is affected significantly by the change in the electron 
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configuration of the (3d) and (4s) shells. As the d state 
population in V was increased, from (3d)^(4s) 2 to (3d)^(4s) i 
the d— bandwidth (e.g. H^, ~ H 12> increased by about 70^. 

The prediction of the nature of the Fermi surface (FS) in 
the vanadium group of metals (i.e. V, l\ib and Ta belonging to 
the group VB) was first maae by Mattheiss [13] himself in 1965. 
Assuming a rigid band model , he predicted the FS of the vanadium 
group of metals by using the results of three different calcula- 
tions on the chromium group (i.e. Cr, Flo and VJ) of metals. Two 
of these were his own calculations for two different potentials 
for W and the third was based on the calculation by Wood [14] 
for iron. It was found by Mattheiss that in contrast to the 
VI B group (i.e. chromium group) metals, the F‘S of the VB 
group (i.e. vanadium group) is rather sensitive to small 
changes in the bands true t ur e . Such a behaviour arises because 
the r 25 , state in the case of the VB metals falls just above 
the Fermi level. The topology of the FS for the vanadium 
group of metals showed two hole sheets:. one corresponding 

to the second band consists of a closeu surface which is 
centred on the point r and is octahedral in character. The 
other sheet, arising out of the third band, consists of a 
jungle-gym ( JG-) structure of hole arms along the <100> direc- 
tions as well as large ellipsoids around the points N. The 
conclusions drawn from the results for Fe, however, showed a 
neck along r IT in the third band; otherwise the three different 



51 


calculations predicted qualitatively similar results. Mattheiss 
also found that the accidental degeneracies in the (100) and 
(110) planes (where the inner closed hole surface contacts the 
outer multiply connected hole surfaces) could be removed by 
including the spin-orbit coupling. 

Snow and Waber [15] calculated the bands tructure and the 

density-of -states curves for the bcc and fee metals in the 

iron-series transition metals using a HSC APW method and full 

Slater exchange (ct=l) . These calculations were performed for 

xx 2, 

two different configurations for the outer electrons, (3d) (4s) 
and (3d) n ’’"' 1 '(4s)' ] ', for all the metals. A comparison with the 
experimental results indicated that the (3d) n+ ^(4s)^ configu- 
ration was generally more favourable. 

Anderson et al. [16] calculated the bands true ture of Y 
using the BBC APW method with the electron configuration as 
(3d)^(4s) for two different values of the exchange parameters 
a=l and a=2/3. They found that as a was increased the d-band- 
width decreased significantly. The p-like point IT-^ was found 
to lie below the Permi level for a=2/3 while it was above E=lp 
for a=l (full exchange). A comparison of these theoretical 
results with those obtained from the photoemission [17] and 
de Haas-van Alphen effect [18] experiments supported the a=l 
exchange. It was also found that their results for a=2/3 
case were qualitatively same as those of Mattheiss [12] 'using 

a full (a=l) exchange and a (3d)^(4s)'*' conf iguration. 

1. f. T. K 1 i; 

CEN ir: -l : ^ 

*«T; 58471. 
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Christensen. [19] studied the effect of changing the 
population of the outer d and s shells on the bandstructure 
of V and Mo. Ihe muffin-tin potentials were constructed 
for the case a-1 by using the He'S self-consistent charge 
densities for the two different configurations and the band- 
structure was calculated by using a non-relativistic NSC AfW 
method. It was found that for 7 the d-bandwidth increased 
dramatically when the 3d population was increased from 
(3d)^(4s)^ to (3d)^(4s)~ whereas similar changes did not 
occur for Mo when the 4d population was changed from (4d)^(5s)^ 
to (4d)^(5s) 1 . The results of Christensen for the (3d)^(4s)* 1 ' 
configuration of V agreed well with those of Mattheiss [12] 
for the same configuration. 

The first SC calculations on V were carried out by Yasui 
et al. [20] using the modified tight binding (MTB) method and 
the OPW method for the d and s valence electrons respectively. 
Using the configurations (3d)^(4s)^ two potentials were 
created, one for a=l while another for a = 0.725. Hie 
results for a=l case were very different from the previous 
NSC calculations while the c: = 0.725 results showed a better 
agreement . 

Self-consistent calculations for V were also reported 
by Papac ons tantopoulos et al. [21] who used the APW method 
and the (3d) (4s) configuration. Calculations were performed 
for the normal lattice constant (a Q ) using a=l and a=2/3 as 
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well as for the reduceu lattice constant (a = 0.95 a- 0 ) with 
<x = 2/3. It was found that their results for a=l a.fter the 
first iteration (NSC) were very close to the results after 
the last iteration (SC) for a = 2/3. Such cancellation of 
the effects of self-consistency ana exchange were also pointed 
out for Ni by Connolly [22] and for Cu by Snow ana Uaber [23]. 
The effect of the different values of the exchange parameter 
was reflectea in the changes in the s,p,d,f and the interstitial 
(i.e. outside the inscribed sphere) components of the charge 
density, lor the SC a = 2/3 case, the d-component was more 
pronounced than for the SC a = 1 case. Hie s-d separation 
and the d and s-p widths were all found to increase in going 
from SC a = 1 to SC a = 2/3. The only important change 
observed in the level ordering was for the state which 
falls below 0^ and 0, as the exchange parameter was reduced 
from a = 1 to a = 2/3. Along ry the third band was found to 
dip below the ii'ermi level when going from a = 1 to a = 2/3, 
thus giving rise to significant changes in the PS. A compa- 
rison with the experimental data from the photoemission [17] 
and dl-IvA [18] studies showed that the SC a = 2/3 calculations 
gave a better agreement compared to the SC a = 1. At the 
same time it was found that their NSC a = 1 calculation 
for the configuration 3d^4s^ gave slightly better agreement 
with experiment than SC a = 2/3. Their calculations for the 
reducea lattice parameter ( a =0.95 a ) showed a significant 
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change in the results for the hand structure and Fermi surface 
of 7. Hattox et al. [24] performed a SC AH'/ hands tructure 
calculation for 7 using a statistical exchange X [5] muffin- 
tin approximation. they assumes a ferromagnetic tcc struc- 
ture as a function of the lattice parameter and calculated 
the cohesive energy, pressure and magnetization for 
several values of the lattice parameter. 

Wakoh and Yamashita [25] have calculated the bandstruc- 
ture of 7 using the state -dependent potentials. Using a 
configuration of (3d)^ °^(4s)' L "' ! ' and the values of exchange 
parameter cc - 2/3, 0.8 and 1.0, atomic charge densities were 
calculated self-consistently to obtain the corresponding 
muffin-tin crystal potential. the modified Uigner-Seitz (KWS) 
potentials were also used with the parameter X = 1.0, 1.4 
and 2.0. The bandstructure was calculated self-consistently 
using the KKR method. It was found that the SC & = 2/3 case 
generally gave better agreement with the experimental results. 
However, the calculated values failed to match the experimental 
results (particularly the dimensions and shape of the FS) 
completely. She calculations were, therefore, done phenomeno- 
logically using a final (i.e. fitted) state-dependent potential 
with different energy dependence for the de and the dy states. 
The results obtained with the final state-dependent potential 
were thus able to reproduce the FS data and the X-ray. form 
factor values obtained experimentally, 'these authors later on 



extended their calculations to compute the directional 
Compton profiles [26] and positron angular correlation 
curves [27] using the same state-dependent potential. We 
shall discuss these results at a later stage. 

Recently Boyer et al 0 [28] have investigated the 
effects of exchange and self --consistency on the bands tructure 
of the group VB metals Y f Eb and fa, fhey have calculated 
the dens ity-of -states , PS areas and effective masses from 
the SC AB; bandstructure results obtained by using Slater's 
[29] X a exchange approximation. In this method the value of 
a is taken to be that for which the HFS total energy of the 
atom is the same as the HP total energy. For V, the value 
of the X parameter previously determined as a = 0.71556 was 
used. A comparison of their results with other theoretical 
and experimental results showed that the SC <x = 2/3 results 
agreed with the experimental results better than their X 
results. It was also found that the bandstructure results 
for V are very sensitive to the exchange approximation used. 

In a recent paper Callaway et al. [50] have reported a 
bandstructure calculation of the paramagnetic 3d transition 
metals (including V) by using the self-consistent 1CG0 
(linear combination of Gaussian orbitals) method. They have 
also calculated the Compton profiles of V and these results 
would be discussed later. 
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Recently Kanhere and Singru [31] Rave reported e calcu- 
lation of the bandstructure and EMU in the 3d metals Y, Cr, 

Re, ITi and Cu by using Hubbard’s fast approximation scheme. 

By using a common method for the construction of the crystal 
potential, the systematics of the EMD and their relation to 
the respective bandstructure were pointed out. 

To summarise, the results of various calculations have 
shown that the bandstructure of V is very sensitive to small 
changes in the crystal potential introduced by (i) the changes 
in the (3d) and (4s) outer electron configuration, (ii) varia- 
tion in the exchange parameter cu Besides these the other fac- 
tor which seemsto affect the bandstructure results is.(.iii) 
whether the calculations were performed in a self-consistent 
(SC) or the non self-consistent (NSC) manner. In Table 2.1 
we have summarised some of these features for those bandstructure 
calculations which show a fair agreement with the experiments 
(and therefore with each other). 

Some of these factors [(i) to (iii) mentioned above] 
affecting the bandstructure results tend to provide a cancel- 
lation effect. These are outlined below . 

(a) Keeping the other factors common, a NSC a = 1.0 
calculation is equivalent to a SC a = 2/3 calculation. This 
behaviour has been pointed out earlier in the literature 
[22,23]. 
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Table 2.1 Summary of the method and parameters used in the 
various bands true ture calculations 


Reference 

Method 

FSC or SC 

Electron 

Conf iguration 

Exchange 

Iiattheiss [12] 

APW 

use 

(5d) 4 (4s) 1 

a=l 

Anderson et al. 

[16] 

Ap'd 

use 

(3d) 3 (4o) 2 

a=l 

Yasui et al. [ 20] 

KTB+OPU 

sc 

( 3d) 4 (4s) 1 

a=0. 725 

Pa paeons tantopoulos 
et al. [21] 

APV 

sc 

(5d)^(4s) 2 

CM 

II 

<3 

u'akoh and Yamashita 
[25] 

. KKR 

sc 

(3d) :> * 9 (4s) 1 “ 1 X a method 
(see text), 
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(b) Using the BSC API method the results obtained by 
Mattheiss [12] for the case a = 1, (3d)^'(4s)^ compared well 

rz o 

those for the case a = 2/3, (3d) (4s) by Anaerson et al. [16] 
Phis suggests that the effect of recoicing the 3d population 
from (3d)^'(4s)^ to (3d)^(4s)^ is compensated by decreasing a 
from 1 to 2/3= 

2 • 3 S urvey of the Fermi Surface studies of Vanadium 

In this section we shall briefly review the experimental 
studies of the Fermi surface of Y. Hie general conclusion of 
the different experiments (which are not many, perhaps because 
of the lack of pure samples of Y) support the FS topology 
suggested by Mattheiss [13,32]= 

The earliest work carried out bj r Alekseevskii and Egorov 
[33] by the magnetoresistance studies was inconclusive, but a 
subsequent work [34] using same techniques reported an 
anisotropy in the IS of Y. Phillips [18] using the technique 
of impulsive fields studied the FS in the {100} crystallo- 
graphic plane. By using magneto thermal oscillation techniques 
Parker and Ealloran [35] extended their previous studies [36] 
and investigated the FS of Y in two planes ClOO} and {llO}. 

Phis is the most extensive FS study of Y and these authors 
[35] have compared their results with other experimental and 
theoretical results . 
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different studies [34-36] agree that the general features 
of the IS of V is as shown in Pig. 2.2. She crosssection of 
the PS in the central {110} and {100} planes are shown in 
Pig. 2.3 as calculated by Mattheiss [12] using energy bands 
for tungsten (fig. 2.3a) ana iron (fig. 2.3b) . 

the closed second-zone hole surface of V centered on r 
is shown in fig. 2.2a. The third zone hole surfa.ce (Pig. 2.2b) 
is quite complicated. It consists of multiply connected 
hole tubes along the <100 > directions (which are known as the 
jungle-gym arms) as well as ls.rge distorted half ellipsoids 
centred at IT. Por the latter type of half ellipsoids, Parker 
and Halloran [35] experimentally determined the semiaxes to 
be 4.632, 4.405 and 3.649 nm“^ along the M-P, PT-T and BF-H 
directions respectiveljr. The measured areas of these half 
ellipsoids [13,35] agree amongst themselves but are higher 
than the API I bands tructure calculations [21]. Similarly, the 
effective masses determined from the dHvA experiments [35] 
gave (nPYm) values which were higher than those calculated 
theoretically [21,28]. 

The model PS of V obtained by Mattheiss [13] using the 
energy bands of Pe predicted that the ellipsoids are connected 
to the jungle-gym ( JG-) surface by narrow necks along the r-N 
directions (Pig. 2.3b). The measurements of Parker and 
Ilalloran [35] indicated that such narrow necks exist along 
the r— PI directions. Por such necks to exist the third band 
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in V along the E direction should not cross the 1'ermi level 
but should lie above it, Bandstructure calculations (for 
the normal lattice constant of V) reported by different 
workers [ 21 , 25 , 28, 30] all show that the third band along 
the £ direction dips below the rermi level and thus preclude 
the existence of such necks. 5a paeons tantopoulos et al„ [21] 
however showed that for the reduced lattice constant 
a = 0.95 and a = 2/3 the ellipsoids connect the JG arms 
and estimated that the contact should occur at an applied 
pressure of 135 kbars , 

In the next section we shall point out that this point 
of disagreement (i.e. existence of necks) between experiment 
and theory can perhaps be resolved by a carefully designed 
measurement of the directional Compton profiles of V. 

2 • 4 Survey of the Compton Profile Studies of Y 

The EMn and the Compton profiles of Y have been studied 
both experimentally and theoretically, ihese studies were 
inspired by the early studies of Weiss and ne Marco [37] who 
measured the X-ray diffraction of Y paired reflections and . 
obtained a e^/t ^ population ratio of about 0.19/0.81 for 
the electron states in V. Using this population ratio Weiss 
[38] calculated the directional Compton profiles of V. This 
calculation indicated an anisotropy in the Compton profiles, 
the difference between J-j_qq( 0) and JpppCO) being about 85S. 



Phillips [39] measured the Compton profiles of single-crystal 
vanadium samples using 17-4 keV Mo-K , X-rays and observed 
that the anisotropy at J(0) was very small. Gamma-ray Compton 
profiles were later measured by ferasaki et al. [40] and 
Paakkari et al. [41] and these studies showed that the aniso- 
tropy was measurable but was smaller than predicted by Weiss 1 , 
model [38]. These experiments showed that the difference 
AJ = [ ~ J-) 1 1 ( 0) ] was about 2 % and the anisotropy 

curves AJ(p ) showed a structure (submaxima and minima) which 
could not be explained by the calculations of Weiss [38] . 

A systematic measurement of the Compton profiles of the poly- 
crystalline samples of V, Cr, Hi, Cu and Zn was reported by 
Manninen and Paakkari [42] but later experimental work [43,44] 
showed that the multiple scattering effects are important. 
Experimental Compton profiles for poly crystalline V corrected 
for the multiple scattering effects showed a close agreement 
[43] with a theoretical Compton profile calculated from the 
PEA model [45] in which only 4s electrons in an outer 

rz p 

electron configuration (3d) (4s) are normalised. 

Wakoh and Tamashita [26] performed a detailed calcula- 
tion of the directional Compton profiles of V using the EKE 
banastructure method and a state-dependent crystal potential 

[25]- A comparison of the theoretical anisotropies AJ(p ) = 

z 

[Jhki(P z ) ~ J h'kn 1 (P z )] obtained by these workers showed a 
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fair agreement with the experimental results of Terasaki et al. 

[40] . fhe spherically averaged Compton profile (due to the 
hand electrons) obtained from these calculations showed a 
better agreement with the measured profile of Paakkari et al. 

[41] rather than that of Phillips [39]- Although Wakoh and 
Yamashita [26] presented the theoretical contributions by the 
first, second and third band to the Compton profiles along 
the <100 >, <110 > and <111> directions and discussed their 
results in terms of the bandstructure and PS of V, they do not 
describe the behaviour of the EMI), P(p) itself. In a later 
work, Wakoh et al. [46] repeated these calculations by using 
the A Ptf method and obtained results which were qualitatively 
similar to their earlier work [26]. 

Recently Kanhere and Singru [31] discussed the systema- 
tics of the EMI), P (p) , In the 3d transition metals V, Cr, Pe, 

Mi and Cu and explained the behaviour of the EMi) in terms of 
the respective bandstructure. However, this analysis was 
restricted only oo the three symmetry directions <100>, <110> 
and <111>. !Phe recent calculation by Callaway et al. [30] 

using the LOGO method concerns itself only with the spherically 
averaged Compton profile from poly crystalline Y. 

One of the main objectives of the present work is to 
gain an insight into the EKE, p(p) , of Y in terms of its 
bandstructure and PS. If the narrow necks along r-N (connecting 
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the JG- arms to the ellipsoids) exist the EMn 5 p(p) s along rilf 
should he radically different and this behaviour should be 
reflected in the directional Compton profiles and their 
anisotropies. ill is aspect would be discussed further in 

Chapters 4 and 5. 
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Chapter 3 


THE APPROXIMATE CALCULATION OF ELECTRON BANDS TRUC TORE 

AND MOMENTUM DENSITY 


3 . 1 Introduc tion 

It may be recalled from Eq . (1.18) that the electron 
momentum density (EMD), p(p), in the independent particle 
model is calculated [1] from 

p(p) = !-> |A..(k,p) | 2 6(p-k“E. ) (3.1) 

where the various symbols are alreadjr explained in §1.4. 
This involves the calculation of the matrix elements 

A.:(k,p) = / e “ ip ‘ r *£ .(r) dr 

J cell J 

4 - *>* 

for each j and (k,p). 

, 

In a typical computation, one has to evaluate the EMD, P(p), 
at several (usually > 500) momenta values of p involving many 

- 4 - -*• - 4 - 

combinations of k and to form p = k + Such computa- 

tions can be very time-consuming if a standard (and more 
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accurate) band structure method like the APW or KER method 
is used to calculate the electron vra . ref unctions i'ft . , the 
banastructure E(k) and the EMb p(p) . It is, therefore, 
useful to employ a rapid but approximate method to compute 
the above quantities. Such a fast approximate scheme for 
transition metals was first proposed by Hubbard [2]. This 
method does not contain any arbitrary parameters and it has 
shown to provide an overall accuracy of about 0.005-0.01 Ry. 
for the bandstructure results. Ill jnar ends [l,3] further 
developed this method for the calculation of the Eft). 
Application of the Hubbard-Mijnarends fast approximate 
method for the calculation of the EH in the transition 
metals like Ee [1], ¥i [4], Cu [1,5], Pd [6] and Ag [6] 
has met with significant success when the theoretical 
resuLts are compared with experiment. In view of these 
successes and the availability of computer programs 35 and 
its suitability and speed for the computer system (IBM 7044/ 
1401) located at this Institute, we have adopted this 
method for the calculation of the EMi), p(p) , in vanadium. 

In the following section we give a brief outline of this 
method, the details of which are described elsewhere [2,7,8]. 


These computer programs were generously made available 
to us by hr. P. E. Mijnarends. We are indebted to him 
for his kind help. 
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3.2 Hub bar d's Fast Approximate Schem e 

In the case of transition metals the electron band- 
structure can be described as hybrids of nearly -free -electron 
(j£FE)”like (or plane -wave -like ) bands with narrow d bands 
of tight-binding nature. Ziman [9] has reviewed various 
approximate schemes using the above concept. One such 
convenient scheme is that due to Hubbard [2] and it is 
based on the Green's function [9] or Horringa [10], Kohn 
and Ro stoker [11] (i.e. KKR) method of bands tructure 
calculations. The Green's function or the ETTR method 
determines the bands tructure through the calculation of the 
energy-dependent phase shifts arising from the scattering 
of electrons from the atomic potential (usually of the 
muffin-tin type) . 

Ziman [12] reformulated the EKR method and obtained 
the determinantal equation 


u 




n 


- e) B + 

' n 

l 

n 1 

r , 

nn 1 



-*■ 

I k ! and 

k = 

k + 

i n 1 

n 


B n' = 0 


(3.2) 


where k n = | k^ | and k n = k + IC, being the reciprocal 
lattice vectors , and e is the Bloch energy eigenvalue for 
the electron in a state k. 


It was shown by Lloyd [13] that B n are the coeffi- 
cients of the pseudo-vavef unction $ _ such that 

ps 
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'P ps (r) = (1 /t) 2 l B n exp (ik n *r) (3.3) 

throughout the atomic polyhedron (of volume t) . In the 
interstitial region outside the inscLioeu spheres, the Bloch 
wavef unction is given by the same expression, i.e. 

%(r) = (l/r) T I B n exp (ik n -r) (3.4) 

n 

.The matrix elements r s are given [14] by 


nn 1 


V l < 21+1 > 4 ntlVP -l {k n> r l ) P l !cos e xm'> 

(3.5) 


where 


-tan 

1 kr i kr f ~( tan 'f n f kr i^ 


(3.6) 


k = V"e, angle between k n and , , r^ is a radius 

in the interval (0,^), r^ being bhe muffin-tin radius j 
is the phase shift of the 1th partial wave scattered from 
the atomic potential while n^ are the spherical Bessel 
and Neumann functions respectively. 


Exce 


pting the case when 3 p(kr-^) or [ j-^(kr^) -( tan p^) 


n-j(kr^)] is close to zero for some value of 1 we stall 


have and therefore I* , to be small. It has been 
suggested by Pendry and Capart [15] that in order to avoid 
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the large values , one 
n-j (Icr-^) if there is a 
be chosen at = r^, 
The case of being 


should choose at the zero of 
zero in (0, r^) , otherwise r^ should 
the radius of uhe inscribed sphere, 
different from r^ Trill occur only for 


1 = 0 . 


It has been shown by Hubbard [2,7] that in the case of 
transition metals there occurs a narrow resonance for tan 
and hence also for for the d bands corresponding to 1=2. 
For the case of a single resonance at k 2 = e Q for 1 = A( X 
being 2 for the d bands of a transition metal) one can 
write [2J 


where 

thus 


r is t 
making 



he width of the resonance 
y£ always small . 


(3.7) 


and y^, is the residue, 


Ziman's determinantal equation (Eq. 3.2) can now be 
written as 



where the matrix elements of the various submatrices are 
written as 


iV mi * ^h ^nn ' 


v u • y<*/ 


V nn- = 


(^) I (21+1) r* y*_ P 1 (cos e nn .) 


(3,9b) 
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h 


nm 

md 


-4s r x ( r A) 2 J^Vj) T xm (k n ) 


a 


m 


4 it r. j 

.4- (t>/ t ) 2~ T i 0; r ) Y* (V ) B 

e 0 "£ w/ ' £ J X V n X' “Xm u ~n ; n 


(3.9c) 


(3.9d) 


with m = -X, -- X + 1, .... x 


In Eq. (3.8) the column matrix ^ describes the 
resonant behaviour while fc- contains the non-resonant part 
of the wavef unctions. 

By following a proceaure suggested by Heine [16 ] , 
Hubbard [2] was able to reduce Eq. (3.8) to an equation of 

~Qb 

lower order. fhis procedure consisted in dividing k n into 
two sets: 

(a) A preferred set ' P ! which consists only of the chosen 
(few in number) fix for which k 2 has values in or near the 

li il 

energy values of our interest. Usually the set P consists 
of a few of the smallest reciprocal lattice vectors. 

(b) A remainder set ’R* which consists of those which 
are not included in the set P. 

By introducing the approximations V ^ = 0 

( R 1 2 

and by neglecting in comparison with (k^ - e), 

Eq. (3.8) can be split up as 
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ana 


( |lf - b H ■ S’--. + fk b 


(3. 1.0b) 


where tlie square matrix ^ is given y 


A 


mm* 


= £ „ 6. 


(4k) 


o mm 


r j (E) b i vb 

n 




k 2 


x W £ n> Y *Xm’^n> 


.(H) 


(3.11) 


the sum £ being carried out over all the vectors of the 
n (R) 

remainder sot R. A prescription to calculate the sum £ 

n 

is given by Hubbard [2] as follows 


I 


(H) 


= 1 


(S) 


4" 


ram* 

r x T k 


+ in - S±(X) ] 


X=2k c r A 


(3.12) 


«(S) 


where p ' is a sum similar to the second term in Eq. (3.11) 
but now over the first few shells of reciprocal-lattice 
vectors except those of the preferred set 'P', Si(X) is 

y 

the sine-integral / (sint/t)dt, %■, is the volume of the 

o 

unit cell in reciprocal space, and (4Kk^/3y) = E, E being 
the total number of reciprocal lattice vectors included in 

is) 

P and in the sum £ ° The deterrainantal equation (Eq. 

3.10a) resembles in form the equation used in model pseudo - 
potential schemes [16,17]. The physical significance of the 
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various subna trices in Eq. (3.10a) is as follows. The sub- 
matrix 'V describes the nearly-free--eleetron (HFE) part of 
the (s-p) conduction band while the submatrix $ (5x5 in size) 
describes the d--bands. The hybrid iz a. t ion between these 
bands is contained in the off-diagonal matrices and Ik 

/t\ n .• 

The submatrices A; and w depend on energy. The energy 

dependence of f] and V can be suitably approximated to 

{!*/{> 

remove the energy dependence of the matrix v in 

[ fl (e) - tl = 0 (3.13) 

which is an equivalent form of (3.10a). The idea is thus 
to reduce the problem to that of diagonalization. 

Hubbard [2] has indicated the procedure for this. In 
outline it consists in setting 


iM’(e) - eJT ~ fiL — e M, 


(13.14) 


such that 11^ anc -'- ’^2 ao no ^ depend on e 9 and writing 


$ „ 

\ C\ 


[J 1 and y 


i A 

rat * 


to obtain 


[ r^i 




0 


(3.15) 


as an equivalent form of Eq. (3.10a). 
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By comparing Eq. (3.15) ana (5.10a) we note that the 
problem has now been put into a secular form which can be 
solved by standard methods of matrix algebra. 


3 • 3 Calcula tion of Electron Wavefunctions and Momentum 
Density ~ - 1 - ■ - 

In the previous section \<re have given a bars outline of 
Hubbard’s fast approximation scheme. In this section we shall 
outline how one obtains electron wavefunction ’/'•^(r) and 
moment-urn density p(p) in the Hubbard -Mi jnarends method 
[ 1 , 5 ]. 


In the case of a muffin-tin type of crystal potential 
the electron wavefunction;, ^^(r), can be expanded in two 
different ways [13] in the two regions. Inside the inscribed 
sphere (r < in) the wavef unction can be expanded in terms 

A f 

of spherical liar monies ( r ) with 1 = (3., m) as follows. 


=■ 


4n 


h) 


1 x 
I 


.1 


R-.(r) 

H^T5~7 Z l Y l (r) 


( 3 . 16 ) 


In the region (r > r^) between tile muffin-tin spheres, 
the electron wave function is expanded in terms of plane 
waves with B n as the coefficients 

*£(?) = [ B n exp[i (k+lt) *r ] (3.17) 

(t) 2 n 


The coefficients Ij- in Eq. (3.16) are determined by applying 
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the continuity condition on at r=r. to obtain 

Lu jL. 


x L = I B n an Vi) ^ 


(3.18) 


One will have to use all 1 -values upto infinity to achieve 
a perfect matching at r = r.. In practical computations 
one truncates the expansion in Eq. (3.16) [and hence in 
ilia. (3.18)] at a finite value of 1 (usually 1 = 3) and 

this gives rise to a slight discontinuity at r=ru. 

In the fast approximation scheme of Hubbard [1-3] the 
expression for A^(k ? p) can be simplified by using Eq. (3.1 6) 

J 

and (3.17) to obtain 


A.j(k,p) 


I B n 6( V p) 


II 



“h 



r . 


i t u » (5) l 1 h ( ^ R i (r)rSa 

1 m 1 v x o 

(3.19) 


In the present calculations for vanadium we have used 
Eq. (3.1) and (3.19) as the basis to compute the EMI), p(p) 
due to the band (3d) (4s) electrons. The MU due to the 
inner or core electron [(Is) (2s) ( 2p) v (3s) ^(3p) ] states 
can be calculated by using Hartree-Eock wavef unctions and 
such results are available in literature [18]. 
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3.4 C omp utationa l Detail s 

The course of the present computations progressed 
along the following stages. 

3.4.1 Calculat ion of th e muffin- tin cry stal potential 
for vanad ium 

The starting point of our calculation was the choice 

of a muffin-tin potential. Using the lattice constant a 

for metallic b.c.c. V as a = 5.74509 a.u. and the radius of 

inscribed sphere r^ = 2.4877 a.u. we have done complete 

calculations for two different muff in- tin potentials for V. 

The first one (called 1/1 hereafter) was taken from the 

calculation of Mijnarends [19]. This potential l/l(r) was 

constructed by Mijnarends by using neutral-free -atom 

3 2 

charge densities for the (3d) (4s) electron configuration 
obtained by a relativistic calculation [19,20]. The 
exchange potential was included with a = 2/3, following 
Caspar, Itohn and oham [21]. Using L'dwuin alpha expansion 
[22] atomic charge densities upto 13 th nearest neighbours 
were superimposed. The potential was not made self-consistent. 
The constant potential between the inscribed spheres was 
obtained as Vl^ = -1.7106 Ry . The potential was construc- 
ted on a Herman-Skillman mesh [23]. The quantity rl/l(r) 
normalised to unity at r = 0 on a coarse mesh is shown in 
Table 3.1. 
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Table 5.1. 

The normalised mull in- -tin 

potential (rVl(r)/2Z) 

SI. 

No. 

r 

in a.u. 

rVl(r)/2Z 

1 

0. 00000 

1.00000 

2 

0.00194 

0.99333 

3 

0. 00339 

0.98659 

4 

0.00778 

0.97297 

5 

0.01167 

0.95940 

6 

0.01945 

0.93283 

7 

0.02724 

0.90743 

8 

0.04280 

0.86063 

9 

0.05837 

0.81859 

10 

0.08950 

0.74545 

11 

0.12063 

0.68374 

12 

0.18289 

0.58433 

13 

0.24516 

0.50756 

14 

0.36968 

0.39711 

15 

0.49421 

0.31822 

16 

0.74326 

0.21103 

17 

0.99231 

0.14012 

18 

1.49042 

0.06129 

19 

1.93853 

0.02303 

20 

2.98474 

0.00132 
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A second type of muffin-tin potential (designated V2) 

was constructed by us following the method of Bhatnagar [24]. 

\Ie employed an electron configuration (3d)^(4s)"^ and full 

Slatei 1 exchange a=l. Our efforts to use a = 2/3 for this 

A A 

calculation cf V2 (configuration 3d r 4s J ~) were not successful 
because of -the computational difficulties in the Hubbard 
scheme. this point would be taken up again in the next 
Chapter when we shall discuss the results of our calculations. 
Bor the case of V2 the neutral free -atom-charge densities 
were calculated employing the wavef unctions of Herman and 
Skillman [23]. I'Owdin's alpha expansion scheme [22] was 
again used to superimpose charge densities upto 13th nearest 
neighbours. Ho attempt was made to make the potential self- 
consistent. The constant interstitial potential was obtained 
as Vq 2 - -1.5836 Ry. Table 3.2 tabulates normalised rV2(r) 
for this second potential. 

Towai-cis the end of our computations we received the 
numerical values of the self-consistent muffin-tin potential 
used by Papaconstantopoulos ot al. [25]. This potential 
(called hereinafter as V3) used the Gas par-Kohn -Sham exchange 
parameter cc = 2/3 and was constructed for the normal lattice 
constant (a = 5.713 S the value used in Ref. 25). 

In Table 3.3 we have summarised some important numerical 
data about the three potentials 171, V2 and V3. As explained 



82 


‘.Cable 3.2 

. normalised muffin-tin potential (rV2(r)/2Z) 

SI. Ho. 

r in a . u. 

(rV2(r)/2Z) 

1 

0.00000 

1.00000 

2 

0.00194 

0.99405 

j 

0.00389 

0.98799 

4 

0.00778 

0. 975 66 

5 

0. 01167 

0.96319 

6 

0.01945 

0.93836 

7 

0.02724 

0.91418 

8 

0.04280 

0.86861 

9 

0.05837 

0.82679 

10 

0.08950 

0.75312 

11 

0.12063 

0.69137 

12 

0.18289 

0.59301 

13 

0.24516 

0.51596 

14 

0.36968 

0.40111 

15 

0.49421 

0.31849 

15 

0.74326 

0.21107 

17 

0.99231 

0.13985 

18 

1.49042 

0.05861 

19 

1.98853 

0. 01813 

20 

2.93474 

-0.00216 



Table 3 . 3 . 


numerical constants for the three potentials 
for V. 


parame ter 

Y 1 

V 2 

t 3 

a, the lattice const, (a.u. 

) 5.74509 

5.74509 

5 .713* 

rpj the muffin-tin radius 
(kept as Jr the nearest 
neighbour distance) (a.u. ) 

2.4877 

2.4877 

2.4738* 

Constant potential between 
the muff in -tin spheres 

(By • ) 

-1.7106 

-1.5896 

-1.5498* 

Ep, the Fermi Energy with 
respect to the zero of 
the muffin-tin potential 
(By . ) 

0.746 

0.774 

0.777** 

e ^ s the mean d band energy 
(the energy at which tanpp 
goes through a resonance 
(By.)) 

0.838 

0.863 

0. 843** 

¥iath of the resonance 
(By..) 

0.212 

0.184 

0.193** 


Taken from Papaconstantapoulos et al. [25] 

Obtained, by us by using the muffin-tin potential of 
Papaconstantapoulos [25] in the approximate method of 
Hubbard [2], 
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in Chapters 1 and 2 the aim of the present work has been to 
examine the sensitivity of the EMu, P (p) and the directional 
Compton profiles P z ) J co the choice of crystal potential, 

electron configuration etc. 

3.4.2 Calcul ati on of Phase Shift s a nd ha dial Wavefunct ions 

After choosing a particular muffin-tin potential for V, 
it was interpolated on an Herman— Skillman mesh [23] using 
£>x = 0.0003125. ‘fhe next step was the calculation of phase 
shifts (see Eq. (3.6)), in the energy range E = 0.0 

to 2.0 Ry j at about 200 energy values situated at an energy 
interval of 0.01 Ry,. Values of R-^(r), the radial part of 
the wavefunction was stored on a magnetic tape for further 
use. 


3.4.3 Calculation of line rgy Hi ge nva lue s ( ) , wave f unc t io ns 
(r.e r% and Xp coefficients) and momentum density 

fhe phase shifts t>^(E) and the radial wave! unctions 
Rl(r) having been determined we now proceeded to calculate 

®n J ®n J arL< ^ eac k selected k -value by using the 

scheme [see Eqs. 3.1 to 5.1?] described in §3.3. 

Within the inscribed sphere the wavefunction expansion 
(Eq. 3.16) was truncated at l max = 3 while in the interstitial 
region 141 plane waves were used (Eq. 3.17) to expand the 
wavefunction. 
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In his original calculation on the 3d transition metals, 

Hubbard [2,3] had usea n = 7 ( n being the number of vectors 

P P 

in the preferred set P) for the bcc structures and n^ = 4 

for the^fcc structures. Later on Kaga [26] compared the 

results of band structure calculations performed on Pd using 

the TiTvR method and Hubbard's method and showed that for Pd 

the number n should be hiaher than 4. The results of 
P 

ICanhere [27] and Harthoorn [8] showed that increasing n 

ir 

from 4 to 14- changed the E(i?) and p(p) values significantly 
for Pd but for Ag the change was within the tolerable limit 
(e.g., aE < 0.01 Ry, which is the accuracy of the Hubbard's 
method) . 

¥e examined the effect of n^ on the E(k) and p(p) 

values obtained for V using the potential 1/1. It was found 

that in going from n = 7 to n = 14 the changes were signi- 

P P 

ficant but increasing n^ from 14 to 19 resulted in less 
serious differences. Some typical results are shown in 
‘Table 3.4. It should be pointed out that the calculations 
for n = 19 were performed at the IBM 370 computer at 
Kernforsch.ungsanla.ge JUlich, West G-ermany [28] which had 
the requisite larger memory. The effect of changing n^ on 
the important energy gaps in. 7 is shown in Table 3.5. The 
results of Tables 3.4 and 3.5 show that changing n from 

ir 

14 to 19 causes energy changes within the accuracy of 
Hubbard's approximate scheme except for the points P where 
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fable 3.4. Effect of the number (n p ) of preferred vectors on 
energy values of V s at some symmetry points. 

(Energy®"-' in Ry . ) 


r 

n 

P 

14 

19 

r l 

0.225 

0.228 

r 25 ! 

0.795 

0.800 

r 25 * 

0.795 

0.800 

r 25 ' 

0.795 

0.800 

J 

r i 2 

0.941 

0.944 

u 

r 

12 

0.941 

0.944 



h 12 

CM 

H 25 ' 

H 25 ' 

H 25 ’ 


H 


14 

0.502 

0.503 

1.053 

1.055 

1.055 

1.462 


IS 

0.505 

o . 503 

1.054 

1.057 

1.057 

1.471 



% 

N 2 


H* 1 
H 1 

) 

E 4 

% 

S 


14 

0.474 

0.670 

0.818 

0.942 

0.971 

1.086 


IS 

0.468 

0.673 

0. 823 

0.944 

0.972 

1.086 



P 4 

P 4 

P 4 

P 3 

p 3 

P 4 

P 


14 

0.659 

0.664 

0 * 656 

0.975 

0.975 

1.400 


19 

0.651 

0 . 654 

0.654 

0-377 

0.977 

1.282 


Potential used vas 1/1. 

Values given with respect to the zero of the muffin-tin 
potential. 
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Table 3.5. hi feet of n. on the energy gaps in 7. 

Ir 




Case 

71 



n v = 14 

n p = 19 

s-d 

r - r 

25’ 1 

0.570 Ry . 

0.572 Ry. 

separation 

0.830 

0.C29 

tt ~ r 

25 f 1 

d 

bandwidth 

H 25' “ H 12 

0.552 

0.554 


r i2 " r 25’ 

0.146 

0.146 

s-p 

bandwidth 

- r i 

0.593 

0.596 


h 15 “ r i 

1.237 

1.243 

Occupied 

banduidths 

% “ r i 

0.521 

0.517 


“ H 12 

o. 244 

0.241 
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Table 3.5. Effect of the number of preferred vectors n.p 
on the WvJ of vanadium. Results for some 
symmetry points it = p. 



n p 

1 

2 

3 

4 5 6 

Total 


14 

0 „ 686 




0.686 

r 

19 

G .578 

- 

- 

- - 

0.678 


14 

0.236 



0.346 

0.238 

H 

19 

0.224 

0.014 

- 

0.344 

0.238 


14 

0.200 


0.385 

0.001 - 

0.200 

N 

19 

0.200 

- 

0.383 

0.002 ~ 

0.200 


14 

0.343 


0.013 

0.088 

0.356 

P 

19 

0.380 

- 

- 

0.077 

0.380 
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fable 3.7. The set of preferred vectors in units of 2 %/& . 

In the calculations only the first 14 were used. 


1 

(0,0,0) 

2 

(0,-l,l) 

3 

( 0,-1, “1 

4 

(0,-2, C) 

5 

(1,"1, 0) 

6 

(~i,~i,o; 

7 

( -1 9 0 y *“"-L , 

8 

(“1,1, C) 

9 

(1,0,1) 

10 

(1,0, -1) 

11 

(-1,0,1) 

12 

(1,1,0) 

13 

(0,1,1) 

14 

(0,1, -1) 

15 

(2,0,0) 

16 

( — 2 , 0 , u ) 

17 

(0,2,0) 

18 

(0,0,2) 

19 

(0,0, -2) 
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the changes are larger. In fable 5.6 ire have shown the 
effect brought in the EMw, p. (p) , due to the individual 

J 

bands by changing n from 14 to 19. Again these changes 
are minor at the points r ? H, IT but are serious at the 
point P. As the memory limitations of the computer (IBM 
7044) at III Kanpur restricted us from using a larger 
preferred set we carried out our calcinations for n^ = 14 - 

The set of the 14 preferred, vectors used in our 
calculations are tabulated in Table 3.7 which also shows the 
vector numbers 15 to 19 (i.e. the set of IS vectors) used 
in testing the effect of n . It should be remembered that 
by using the set of only 14 vectors we have calculated the 
EMn within a certain accuracy. However, the inaccuracies 
(or uncertainties) brought by this limited set is within 
the limitations of the Hubbard -Hi jnarends method. 

The results of our calculations for the cases VI, 

72 and V3 are discussed in the following chapters. ¥e have 
used about 80 h-points in the l/48tli part of the Brillouin 
zone (Pigs. 3.1 and 5.2) and by choosing suitable Kw" 
vectors we have calculated p(p) for about 700 -points in 
the (100) and (110) planes. 





points in the (100) and (110) planes involved in the 
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Chapter 4 


ELEC (EROS' MOMENTUM -DISTRIBUTION IN VANAuIUM : 
RESULTS ANN ilSCUSSIOi: 


4*1 Introduction 

In this Chapter we shall present and discuss the 
results of our calculations oi the bandstructure , Fermi 
surface (FS) and electron momentum distribution (EMI)) in 
metallic V. The starting point of the calculation was a 
choice of the muffin-tin crystal potential for metallic V. 

As explained in §3.4«1? we chose three types of potentials, 
the basic parameters of which are summarised in Table 4.1. 
Perhaps it might have been more appropriate to choose a= 2/3 
while constructing the potential V2 so as to examine the 
effect of changing the electron configuration alone. Although 
such a crystal potential for the configuration (3d)^(4s)"** and 
a = 2/3 was constructed by us, we met with serious computa- 
tional difficulties in calculating the bandstructure using 
such a potential. These difficulties arise out of the 
following considerations. 

As pointed out in §2.2 the width of the d band in V 

3 2 

is increased when the configuration is changed from (3d) (4s) 
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'Table 4-1- nifferent muffin-tin crystal potentials 

(for metallic 7) used in the present work. 


Our 

label 

Reference 

Electron 

configuration 

Exchange 

parameter 

Lattice 
constant a 
(a.u. ) 

71 

Mijnarends [1] 

(3d) 3 (4s) 2 

a = 2/3 

5-74509 

72 

-T7* 

Present work” 

(3d) 4 (4s) 1 

/V — 1 

u. — JL 

5.74509 

73 

Pa pa c o ns tant o - 
poulos [2]9 

(3d) 3 (4s) 2 

a = 2/3 

5-713 


35 The method of construction of this potential has been 
given in §3.4.1 and. is described by Bhatnagar [3]. 

^ This crystal potential has been constructed in a self- 
consistent manner. 
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to (3d)^(4s)^. Similarly the width of the d hand is 
increased by decreasing a from 1 to 2/3 keeping the configu- 
ration same. These considerations suggest that in going 
from the case VI (where a = 2/3? and (3d; (4s) ) to the case 
where a = 2/3 and the configuration is (Sb)‘ T (4s)^ the d band 
would become much broader. The fast approximate scheme of 
Hubbard [4] is based on the main assumption that the d bands 
are narrow. In the case of the 3d metals the d bandwidth for 
Cu is 0.27 Ry. while for V it is 0/0.55 Ry. [5]. Thus the 
case of V falls just near the limit of Hubbard's approximation 
scheme. If the d band in V, for some crystal potential s 
becomes much broader than 0.55 Ry. the d band becomes more 
delocalised and the main assumptions underlying the scheme 
of Eq. (3.2) -(3.19) break down, giving rise to computational 
difficulties. In fact our attempts to use Hubbard's fast 
approximate scheme to calculate the bandstructure of Hb for 
the case a = 2/3, (4a) (5s) (which is known to give d band- 
with'V'0.75 Ry. ) also met with similar uilf iculties . 

To surmount these computational difficulties we used 
a value of a = 1 in constructing the potential VI. In going 
from VI to V2 the increase in the d population and the 
increase in a have opposing effect on the d bandwidth and 
in fact our results [Table 4.2] indicated that the d 
bandwidth for the case V2 was narrower (by »v 0.04 Ry*} when 
compared to that for VI. 
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The present results for the EM L> and Compton profiles 
(to "be discussed later) obtained from the potentials VI and 
V2 show similarities and suggest that experimentally one may 
not be able to distinguish between the cases VI and V2. 
Recently, a comparison of the experimental Compton profile 
data for some 3d metals has been made with the REA model 
calculations [6] to conclude that a particular electron 
configuration [o.g. , li:(3d) (4s) , 7: (3d) (4s) ] is supported 
by the Compton profile data. She results of our calculation 
would show that for a more realistic comparison with the 
experiment, a band structure calculation (involving a choice 
of the electron configuration as well as the exchange 
parameter) is necessary. It will be shown that although 
crystal potentials can be different, the electron configura- 
tion and the exchange parameter a can combine to give 
similar results for the E ME. Besides this effect, the 
sensitivity of the EHJJ in V to the bandstructure and FS wi ll 
be discussed in detail. 

4 . 2 Energy Bands and Bandwjdths 

The bandstructure diagrams obtained by us along the 
principal symmetry directions of V using the crystal poten- 
tials VI, V2 and Y3 are shown in Figs. 4. 1-4. 6. There is 
a close similarity in the results for the cases 1/1 and V2 
while the results for V3 are also similar except along the 



Case VI 









principal aynuiie 






principal 






m* * -nt'i 





Case V 3 



of vanadium along the prinoip 
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Table 4.2. Comparison of energy bandwiaths (in Hydbergs) for V 


Method and 
reference 

s-a 

separation 

d bandwidth 

s-p bandwidtl 

r 25 »“ 

r l H 25 i_r l 

H 25’- H 12 

r -r 

12 25’ 

h-h 

H 15 r 3 

APW 

Papac ons tan to - 
poulos et al. 
[2] 

0.543 

0.791 

0.535 

0.140 

0.583 

1.220 

rooi 

w'akoh and 
Yamashita [7] 

0.521 

*— 

0.437 

0.157 

0.580 

— 

Present work 
case VI s 

0.570 

0.830 

0.552 

0.146 

0.593 

1.237 

Present work 
case Y2 

0.554 

0.793 

0.515 

0.134 

0.585 

1.250 

Present work 
case V3 

0.514 

0.768 

0.531 

0.142 

0.608 

1.277 

Present work 
case V4 SS 

0.446 

0.653 

0.438 

0.120 

0.565 

1.228 


These results are different from those reported by Mijnarends 
[1] because the value of np (number of preferred vectors) in 
the two calculations are not same. See §3.4.5 for details. 

JBi (3<l) 3 (4s) 2 i a = 1. 
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Sable 


a) r 


b) H 


c) IT 


d) P 


4.3. Enex-gy values for the various symmetry points 


(ener 

■gies given in 

Rydbergs 

with respect 

to the 

zero 

of the muffin- 

■tin poten 

tial) . 


Preference 

r l 

r 25 ' 

r 12 

B P 

71* 

0.225 

0.795 

0.941 

0.746 

72* 

0.263 

0.817 

0.951 

0.774 

73* 

0.288 

0.802 

0.944 

0.777 

Ref .[2] 

0.269 



- 


Reference 

H 12 

h 25 , 

%5 

71* 

0.502 

1.055 

1.462 

72 s 

0.541 

1.056 

1.513 

73* 

0.525 

1.056 

1.563 

Ref .[2] 

0.525 


— 


Reference 

2I 1 

JT IT 

U 2 L 1 

% 

*4 


71* 

0.474 

0.670 0.818 

0.942 

0.971 

1.086 

72* 

0.500 

0.702 0.348 

0.950 

0.979 

1.085 

73* 

0.493 

0.684 0.896 

0.944 

0.973 

1.087 

Ref .[2] 

0.523 

0.692 - 

-- 

- 

- 

Reference 

? 4 

P 3 




71* 

0.664 

0.975 




72* 

0.680 

0.982 




73* 

0.680 

0.976 




Ref. [2] 

0.680 

— 





Present work 
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l direction and their general features compare well with, the 
band structure of V reported by papaconstantopoulos [2] using 
the APV. r method and by Wakoh ana Yamashita [7] using the KER 
method. A aetailed comparison between various results is 
made in terms of the bandwidths in Table 4.2 and in terms of 
the 25(k) values in fable 4.5. Comparing the results for the 
cases 71 ana 72 we note that the potential 71 produces larger 
d bandwidth and s-d separation. Potential 73 resulted in 
smaller s-d separation. Je have also calculated the d band- 
width for the fourth type of crystal potential 74 in volving 
the case (3d)"'(4s) , a = 1 (results shown in -Table 4.2). These 

results indicated that the increase in the d bandwidth in 
going from (3d) 5 (4s) 2 ? a = 1 to (3d) 4 (As) 1 , a = 1 is smaller 
than in going from (3d)^(4s) 2 y a = 1 to (3d) J (4s) 2 , a = 2/3. 

It is thus seen that the bandstructure of 7 responds more 
sensitively to the changes in the exchange parameter a than 
to the changes in the 3d population. We found that the 
d bandwidth increases in the following order-; (i) (3d)^(4s) 2 , 
a = l s (ii) (3d) 4 (4s) 1 ; a = 1, (iii) (3d) :j (4s) 2 , a = 2/3 and 
(iv) (3d) 4 (4s) 1 s a = 2/3. In the last case the d bandwidth 
goes beyond the limits of the applicability of the Hubbard's 
approximation scheme. 

Comparing the results (Table 4.2) for the same potential 
(73) obtained by the AP¥ method (2) and by the present calcu- 
lations we find that although the d bandwidth is about the same s 

4 
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our results show a smaller s-d separation ana larger s-p 
bandwidth. '.These differences are perhaps arising out of the 
aifferent accuracies of the A method and Hubbard's fast 
approximation scheme. 

4 * ^ Fermi Surfa ce, of Vanadium 

In view of the limited accuracy (0.005—0.01 Hy . ) of 
our bandstructure calculation we adopted an approximate 
method due to Burdick [8] to determine the Fermi energy. Bp. 
Using a three-dimensional mesh of 56 k-points in l/48th part 
of the Brillouin zone we calculated the E(k) values for the 
lowest 6 bands. These energy values were arranged in an 
array of ascending order, by giving each k an appropriate 
weight arising from the number of like vectors H(k) for 
each S. Starting from the bottom of the array the energy 
values E(k) of the [TT(k) x-g-xn^jth and the next state (where 
n fe = number of band electrons, n^ = 5 for Y) were found. 

'The average of these two values was taken as tiie Fermi 
energy E-p and is shown in Table 4.5. The dens ity-of-states 
histogram (using AE = 0.02 By. ) for the case VI is shown 
in Fig. 4.7 and it shows the same general features as that 
obtained by a detailed calculation [9]. Using these values 
of Ep and our bandstructure results we obtained the Fermi 
surface of Y graphically. The central crosssection of the 
FS in the (100) and (110) planes of Y for tire cages 1/1, V2 
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ana V3 is shown in Figs. 4.8, 4.9 and 4.10 respectively. 
Except for some minor differences the PS results for the 
cases VI, V2 ana V3 agree with each other and display the 
familiar features (i.e. r -centered octahedral hole surface, 
jungle-gym hole arms and IT-centered hole ellipsoids as 
described in Chapter 2) common to the FS of the group VB 
metals. When compared to the results obtained by Mattheiss 
[10] using the bands for W and by Wakoh and Yamaghita [7], 
the present results show similar topology except for the 
accidental degeneracy in the (100) plane (between the second 
and third band) which is not observed by us. In this respect 
our results agree better with the F5 topology obtained by 
Papaconstantopoulos et al. [2]. 

We have not calculated the FS areas of various orbits 
etc. as we were not primariljr concerned with the de Haas— van 
Alphen data. However, a rough estimate of the size and the 
shape of the ellipsoid centered on IT can be obtained from the 
ratio of the semiaxes along the 1T-P, IT- rand H-H directions. 
We obtained the ratios N-r/H-P = 0.92 and W-r/H-H = 1,34 for 
the case VI and IT-r/H-P = 0.87 and H-r/H-H = 1.33 for the 
case V2. For V3, the corresponding values were 0.90 and 
1.40 respectively. These values may be compared with the 
values H-r/H~P = 0.946 and E-T/H-H = 1.21 obtained experi- 
mentally by Parker and Halloran [11]. In view of the 
approximate nature of our method, these differences are not 


so serious. 
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It may be pointed out that in all the three cases VI, 
V2 and V3, we observe the higher 2-^ band dipping below the 
Fermi level (Figs . 4.1, 4.3 and 4.5) by an amount varying 
between 0.010-0.020 liy. iChe present x’esults for the three 
bands true tures , therefore, do not support the possibility of 
narrow necks (along r~lT) connecting the jungle-gym hole arms 
to the ellipsoids (see §2.3). 

4 . 4 Electron Momentum distribution in Vanadium 
4-4.1 Introduction 

We shall now discuss the results of the electron 
momentum distributions (EMn) obtained by us for the crystal 
potentials VI, V2 and V3, in terms of their bands tructure. 
Such a discussion becomes more meaningful when it is carried 
out on the basis of the symmetry direction involved and 
towards this aim we shall present our results according to 
the symmetry direction (i.e. <100>etc.) involved. 

As discussed in §1.4.3 the total EML> 

. occ + 

p(p) = l p (p) ( 4 . 1 ) 

j J 

(where j is the band index) at a point p = $ + K3 is deter- 
mined by |A^(k,p)| [see Eq. (1.18)]. She magnitude of 
|Aj(k,p)( is, in turn, determined by the nature of the band 
wavefunction j(**)» through Eq. (1,19), as well as by the 
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reciprocal lattice vector involved, in constructing the 

4 

momentum vector p = k + E. . A third important effect on the 

|A-(k,p ) \ was not discussed in Chapter 1 and it arises from 

the group theoretical selection rule pointed out by Mijnarends 

■> 

[12,13] . To understand, this rule one writes A-(k,p) as 

J 


A i (k,p) = N J 
J i 


o( p-k-K. ) / e “ip*z* ^ (r) dr 

1 cell J 

(4.2) 


The symmetry properties of A.(k,p) are determined by the 

J 

behaviour of A . under the operations R of (^(k) , the group 
of k [14]. The function ^ . is the basis function for the 

irreducible representation H of ^T Q - The integral in 
Eq. (1.19) will then be non-vanishing only if the term 
exp (-ip-r) contains a part which transforms according to 
the same representation. To test this we define a projection 
operator 


g> 3 = (i,/h) I x 3 (R) a p. 


R 


R 


(4.3) 


where 1. denotes the dimension of r J , h is the order of 
0 , 

£^ 0 (k), X^(R) indicates the character of the operation R for 
this representation and is the operator corresponding to 
R« The projection operator (P^ , by operating on a function 
will project out of ^ that part which belongs to . 
If it so happens that (p *■* exp( -ip* r) = 0, the exponential 
term does not contain such a part and A-(k,p) = 0 for this 
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particular p« If however & ^ 


exp(-ip-r) ^ 0 the symmetry 


considerations allow A.(k,p) to be non -vanishing . Applica- 

J 


tion of this selection rule is very important while analyzing 


the effects of bands tructure on the BMP due to the individual 


bands 


The above selection rule when applied to the bcc 


Q 

structure (cubic space group Or) gives the following interest- 


ing results [13]. tor a p-direction along the <100> direction 


(if p = k in the first zone or p is parallel to k in a higher 


zone) only the bands of A^ symmetry will give a non-vanishing 


contribution (during its occupied part) to A.(k,p) while the 

J 


other A bands (e.g. A 2 , A ^ , etc.) will give A.(k,p) which is 
iaentically zero. Along the <110> direction only the ^ 
bands and along the < 111 > direction only the band (between 
Tand P) and band (between P and H) may contribute to 


A-(k,p). A more detailed discussion of the effects of this 
3 


selection rule is given by Harthoorn and Kijnarends [13] and 


we have summarised in Table 4*4 some of their results which 


are relevant for our discussion. 


To summarise, therefore, the partial contribution to 


the momentum density at momenta p, due to an electron belong- 


ing to the jth band is determined by the following considera- 


tions : 


♦ 4 * 

(i) The relation p = k + E- 


( ii) Whether the jth band is occupied or unoccupied at k. 
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Table 4.4. Some 51 energy bands allowed by the symmetry 
selection rule to contribute to the Slid at 
momenta p for the bcc lattice. Based on Ref. 

[13]. 


o 

■y 


Observable bands 


z 


0 0 5 

0 € 5 

4 £ 5 

0 5 1 


\ 5 ^1 51 

p l ? P 4 S 


E 12> H 15 


*i 


S] C 5 H 9 , , Hoc: » 


I 1 


15 5 2 * 5 25 

1 


(rm) + 


5 rj n (rHP) +s (r srp) + 

£ /o, r) ^ 0 


“list not exhaustive. See Ref. [13] for a more exhaustive 
tabulation. 
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(iii) Mature of the wave function ^ t ■ (r), 

^ j J 

(iv) Whether the selection rule allows the jth band, (having 

+ + 

a particular symmetry) to contribute to A_-(k,p) or 

J 

not . 

vie shall apply the above considerations while analyzing the 
results of the EiiD . 

4.4.2 Ei'tu along the < 100 > direction 

The EMD, p(p) , along the <100> direction in V is shown 
in Figs. 4.11, 4.12 and 4.13 for the potentials VI, V2 and V3 
respectively. 'The relevant bands tructure (along the a direc- 
tion) is shown in the inset of each figure. Same units and 
scales are used in the p(p) and E(k) plots throughout. It is 
observed that VI, ¥2 and V3 which result in similar band- 
structure also yield similar results for the Mb. As pointed 
out in the previous section only the A^ band is allowed by 
the selection rule to contribute to the EMi) [ through A . (k, p) ] . 

J 

In all the three cases VI, V2 and V3 the higher A^ band 
remains unoccupied throughout as one goes from r to H 
along A. She total EMU, therefore arises only from the 
lowest A-^ band. Since this band never crosses the Fermi 
level, there are no sharp discontinuities in the EMM curve 
which shows a smooth momentum dependence. It should be noted 
that there is a crossing of the second band (from the bottom) 




ari 








The Ml) alon;' the <100> direction for vanadium (case 12 ). 
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with the Fermi level but the symmetry of this band being 
it is forbidden by the selection rule to contribute to the 
EMu . 

Along the < 100 > direction the F8 of V is known to 
exhibit (Figs. 4.8-4.10) interesting xeatures, that is, the 
r~centerea octahedral hole surfaces and the jungle-gym ( JG-) 
hole tubes (arms). It may appear surprising that the EMI) 
along the < 100 > direction (Figs. 4.11-4.13) is relatively 
smooth and aoes not reflect these interesting features of the 
FS. To appreciate this one should remember that these two 
sheets of the FS in V arise from the second band (symmetry A ) 
and the third band (symmetry A,_). fhe symmetry character of 
both the second and third band is such that they are forbidden 
by the selection rule to contribute to the EMI). On the top of 
this the third band is totally -unoccupied. In other words, 
it is only the filled electron surfaces (and not the empty hole 
surfaces) which can contribute to the EEL) provided that they 
are allowed to do so by the symmetry selection rule. 

It would be interesting to further examine these effects 

of the symmetry selection rule by looking at the EMb, p(p), 

along a p-direction which is removed from the < 100 > direction 

but still lies in the (100) plane. One such p-direction 

marked by A in Fig. 4.14 is such that p = (p . p . p ) = 

x y z 

Fhe EMI), in the first Brillouin zone (BZ) along this p-direction 
is shown in Fig. 4.15 along with the relevant bands tructure. 




EMD P(p) 
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It will "be observed that the second band has now lost its 
higher symmetry ( A 2 ) and has assumed a lower symmetry which 
is now (+) . Application of the symmetry selection rule 
(Table 4.4) allows the (+) bands from the THE plane to contri- 
bute along the <£ , r) s 0 > direction. This effect is seen 
from the partial EMn contributed bjr the second band (during 
its occupied part) to the total Slid. 

4.4.3 EMI along the < 110_> directio n 

The results for the EMu along the symmetry direction 
<110 > for the cases ¥1, V2 and V3 are shown in Pigs. (4.16)- 
( 4 . 18) along with the corres ponding bandatructure . When 
compared to the < 100 > directions these results show a more 
interesting behaviour. Along the <110 > direction, the sym- 
metry selection rule allows only the E^ bands to contribute 
to the EMi>. Besides the lower [marked by (a) in the 
bands t rue ture of Eigs. (4.16) -( 4 . 18) ] the next higher E^ band 
[marked by (b) in I'igs . ( 4 . 16) --(4.13) ] can also contribute 
while the band (second band from the bottom) is forbidden 
by the selection rule. In our banastruc ture results the 
higher E^ (or the 1 b * ) band dips below the Eermi level and 
thus remains unoccupied except in the central region along ‘ 
r~E. It is only in this central occupied region that band 
( b) contributes to the total EMh (see the curves of the 
partial EMU markeu by (a) ana (b)) and this explains the 



<11 
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'spikes' and the sharp discontinuities seen in Figs. 4.16- 
4.18. The thinly drawn curves (a) and (b) in these figures 
denote the partial BMu, p.(p), while the thickly drawn curve 
shows the total ERL), p(p) = I P i (p). The behaviour of the 

j J 

partial EMu curves provide a good insight into the combina- 
tion of effects mentioned at the end of §4.4.1. It will be 
noticed that the dipping of the band below the Fermi level 
is shown as the discontinuities in the total HiL) curve not 
only in the first BZ (p^.35 a.u., 0.51 a.u.) but also in 
the higher (p*vl.06 a.u., 0.90 a.u.) sones through the 
Umklapp images of the FS (Fig. 4.16). The width of the spikes 
in Figs. 4.16-4.18 is determined by the amount of ’dipping' 
shown by the band while the height of the spikes will be 
decided by the momentum eigenfunction ofV' .(r) through 

, . 11 S’ J 

-*• p 

jA-(k,p) j . In terms of the FS topology, as one goes from r 

J 

to IT, the first crossing* of the ^ band with the Fermi level 
determines the dimensions of the JG hole-tube arms while the 
second crossing determines the semiaxis of the hole ellipsoid 
centered at IT. I± tho model of the Fb in which the JG hole 
arms and the hole ellipsoids are connected by a narrow neck 
(narrow neck model) is valia, the corresponding bandstruc ture 
could be somewhat similar to that shown in the inset in 
Fig. 4.19. We have calculated the Etilb (along < 3.10 > direction) 
for such a model by lowering the Fermi level by about 0.05 Ry. 
from its position for the case VI. These results are- shown 




■•on for Fermi level dipping below the 
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in Fig. 4.19. A comparison of the results in Figs. 4-19 
and 4.16 show that the EKL> along the <110 > direction changes 
drastically if the narrow neck model is accepted. 


Mijnarends [15,16] ana I'i nolle r [17] have proposed 
inversion methods to reconstruct the H-lr, p(p), by inverting 
(or unfolding) a set of experimentally measured Compton 
profiles ( P z ) ^' or several directions (hkl) through the 

equation 


J(P 2 ) = // P(P) d P X dP y 


( 4 . 4 ) 


Such an unfolding or inversion procedure has not yet 
been reported for the airectional Compton profiles measured 
from V single crystals. he feel that such an unfolding of 
the p(p) from the experimental data for V should be under- 
taken so as to obtain, for example, the p(p) along the <110 > 
direction. The results of our calculations (Figs. 4.16 and 
4.19) would then provide a basis to indicate whether the 
Compton profile data supports a PS model with or without the 
neck. Further in case the experimental results' rule out the 
existence of these necks, the positions (p-values) of the 
discontinuities (measured from r or FT) observed in the 
unfolded p^^^(p) curves can provide some information about 
the dimensions of the JG- hole arms and the N-centered hole 
ellipsoids . 
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ive have also examined the effect of the symmetry selec- 
tion rule on the second hand from the bottom. Along the < 110 > 
direction this band has a 2^ character and is forbidden to 
contribute (Table 4.4) as borne out by Tig. 4.16. If the 
p-direction is moved away from the <110 > direction but still 
limited in the (110) plane the 2^ and ^ bands assume a 
lower symmetry (+) ana (-) respectively. The selection rule 
(Table 4.4) allows only the (+) bands in the (2 Hr!) plane to 
contribute for p = ( K , r], C) direction. Thus the second 
band gives no contribution to the Si even along the p = 

( £ j r), 0) direction (marked by B in iig. 4.14) as shown by 
the results of our calculation (Pig. 4.20). 

4-4.4 EMu along the < 111 > direction 

The plots of the EMLu along the <111> direction in V 
are presented in Pigs. 4.21-4.23 for the crystal jootentials Yl, 
V2 and Y3 respectively. In a bcc structure the direction <111> 
proceeds along r-A-P-F-H-P-P-A-r in the succeeding zones and 
we have therefore shown the bands t rue ture along the A and P 
directions in the inset of Pigs. 4.21-4.23. The symmetry 
selection rule (Table 4.4) allows only the A and F 1 bands 
to contribute to the EMi along the <111> direction and this 
effect is observed in our results. The sharp discontinuities 
in the EMl curves observed at p 0*1.5 a.u. and p 0*2. 25 a.u. 
arise from the crossing of the F^ band with the Permi level 
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the <111> direction (case VI) . 
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and hence it is seen in the second and third BZ. In the 
first BZ no such discontinuity is observed because the 
band is occupied throughout this zone. The position (i.e. 
p-values) of the discontinuities (at p^l.5 a.u. and 2.25 a.u.) 
can provide us with the oilmens ions oi the JG hole arms along 
the H-P direction. As pointed out in the previous section, 
the experimental Compton profile data, if unfolded to obtain 
p 11i(p) 9 can kelp as to determine these dimensions of the PS. . 
The P(p) distributions for the three cases 71, V2 and 73 
(Pigs. 4.21-4.23) show similar features. 

4.4.5 E MU along some other directions in the (1 0 0) and 
TllO) planes "" 

In order to get a better insight into the anisotropic 
nature of the EMU in V we have calculated the p(p) values for 
the case 71 along nine other (besides the <100>,<110> and 
<111>) p-directions lying in the (100) and (110) planes. 

These twelve p-directions are shoxfn in Pig. 4.24 where they 
are marked as 1,2... 12. The components (p , p , p ) of each 

X V Zr 

p vector are shown in Table 4-5 where the (p^, p , p ) refer 

x y z 

to the intersection of the p-airection with the first zone 
face. Each of the p-directions is formed by joining the 
point r to the point- (p^, p y , p^) . 

The plots of the EMu for the directions 1, 6 and 11 
(i.e., the principal symmetry directions) are already shown 




Plane 
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e 4.5- Specification of the components of the twelve 
p directions shown in Pig. 4.24. the values 
(P z »Py»P 2 ) denote a point on the zone face. 


a^J^l of the p 
ction 

4.24) 

P 

in units of 2tc/a 


p y 

P 2 

1 

0.5 

0.5 

0.0 

2 

0.6 

0.4 

0.0 

3 

0.7 

0.3 

0.0 

4 

0.8 

0.2 

0.0 

5 

0.9 

H 

9 

o 

0.0 

6 

1.0 

o 

o 

0.0 

7 

0.9 

0.1 

0.1 

8 

0.8 

0.2 

0.2 

o 

J 

0.7 

0.3 

0.3 

Clo 

0.6 

C.4 

0.4 

Cli 

0.5 

0.5 

0.5 

^L2 

0.5 

0.5 

0.2 
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in Pigs. 4.10, 4-16 and 4.21. Similar plots (however without 
the hancis tructure diagrams) for the other nine p-directions 
(for the case Vl) are presentea in Pigs. 4.25-4.33. All 
these results when examined together illustrate how the 
p (p) distribution changes with the direction of p in the 
(100) and (110) plane ana bring out the anisotropic nature 
of the BMjj. These anisotropies can he understood in terms 
of the bandstructure of V, that is, in terms of two effects, 
(i) the Permi surface topology and (ii) the symmetry selec- 
tion rule effect which allows only the bands of the (+) 
symmetry to contribute along the (g , r) , Oj S ^ r) £ 0) and 
(c» vi, 'I*, { / ti ^ 0) directions. To illustrate this 
effect we have shown the partial EMjj curves for the two (+) 
bands involved in Pigs. 4.25 to 4.33. The vertical arrow 
marked by BZ in these figures indicate the boundary of the 
first BZ. These results can be understood and analyzed in 
terms of the above two (PS and selection rule) effects in a 
manner similar to § 4 . 4.2 ana 4 . 4 . 5 . 

4.4.6 Isodenses of the BUI) in the (100) and (110) planes 

in V 

In order to obtain a better understanding of the EMD, 

■*>* 

p(p), it is useful to plot the results in terms of constant 
Mi) contours (or isodenses showing, the curves p(p) = constant) 
in the (100) ana (110) planes. Besides bringing out the 
anisotropy of the EMU, such plots of the isodenses offer the 










152 


following advantages: (i) She effects of the FS topology «, 
particularly on the high moment-urn components of the EMb is 
more clearly seen. (ii) Shey provide a wider basis to compare 
the results of the unfolding of experimental Compton profiles 
as these results are usually expressed in terms of such 
isodenses [15] • 

Us in g the present results for the EiiD (for the case VI 
and V2) along the twelve p-directions (Sable 4.5) isodenses 
for the total EMb were constructed and are shown in Figs. 4.34 
and 4.35. Shese isodenses provide a wealth of information and 
it would be interesting, in future, to compare the results 
of Figs. 4.34-4.35 (for the total EMb) with those obtained by 
the unfolding procedure . 

It is observed from F'igs. 4.34 and 4.35 that the iso- 
denses of the total Situ obtained from the two crystal poten- 
tials VI and V2 are very similar except for some very minor 
differences. In other words, the isodenses of total EMU 
obtained from the unfolding of the experimental data will not 
be able to distinguish between the cases VI and V2. This 
point would be further discussed in the next section wherein 
we summarise our results . 

From the point of view of the bands true ture, it is 
equally interesting to examine the isodenses of the partial 
EMx) contributed by the individual bands. Ve have plotted such 
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isodenses for the first three (numbered, according to increas- 
ing energy, i.e. lowest is termed first etc.) bands of 
vanadium (case VI ) in the (100) and (110) planes in Figs. 4.5 
4.57 and 4.58. In these figures the (-:-) and (-) signs indica 
the symmetry of the bands in that region and the thick curves 
separating the ( + ) and (-) regions mark the lines where the 
adjacent bands touch or cross. ihere is a transfer of EMJJ 
across these boundaries as tie symmetry of the bands (+/-) 
changes. As seen above only the bands with (+) symmetry 
contribute to the Bi'-ii). In Figs. 4.57 and 4.58 the Fermi 
surface is also markeu, the shaded region denoting the unoc- 
cupied region. 

The only isodenses of the El'-hJ for V reported in the 

$ 

literature are those by Ivakoh et al. [18] for the first band 
in the (100) plane. Our results (Pig. 4.56) show a good 
agreement with their results. These isodenses show consider- 
able anisotropy even for the lowest band in the (110) plane. 
Since this band has almost no crossings with other bands and 
lies wholly in the occupied region below the Fermi surface, 
the anisotropy observed in this case arises mainly due to 
the intrinsic nature of the band itself, i.e., the nature of 
the wavef unctions which describe the band. 

Figs. 4.37 and 4.38 bring out the validity of the 
selection rule, dependence on the FS : topology and the aniso- 
tropy of the partial Mi). They also show that the sheet of tl 
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second band FS (r -centered octahedral hole surface) is only 
partially observable in the (100) and (110) plane by Compton 
scattering. 'The IT-centered hole ellipsoids (third band FS 
sheet) are observable (i'ig. 4.58) only in some places in the 


(100) and (110) planes but with a weak intensity. Same 
appears to be the case with the JG-arus ol the FS. Of course, 
these FS sheets will be more prominently observable outside 


the symmetry planes 
determine the shape 


UOU) and (110) anc. would therefore 
of directional Compton profiles. 


The results of i'ig. 4.57 and 4.58 show the existence 
of nodal lines along the < 100 > and < 111 > directions. As 
explained bp Harthoo'rn and Mi jna rends [13] this is the 


effect of the symmetry selection rule. 


Su mmary 

It is well known that the bands tructure of vanadium is 
sensitive to small changes in the starting potential used 
for the calculation [2,9,19]. To study the effect of 
changing the potential on the bandstructure and EMU, we 
constructed various potentials with different 3d, 4s 
population and the exchange parameters a. However, band- 
structure calculations could not be completed for all of 
them because of the limitations of the Hubbard's 
fast approximation scheme. 


In particular, it was found 
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that the computational method worked s access fully for only 
those potentials for which the d bandwidth was below a 
certain value , namely , vO„55 Ry. 

Prom the results of the potentials studied the follow- 
ing facts seem to emerge; 

(1) There is a variation in the bands tructure when the 3d, 

4s population is change a, the d bandwidth increasing 
with the increase in d population. 

(2) There is a variation when the exchange parameter a is 

changed, the d bandwidth increasing when a is reduced. 

This variation was found to be more significant than in ( 

That is, when a was reduced from 1 to 2/3, the increase 

in the d bandwidth was more, than for the case where 

3 2 4 1 

the d population was increased from 3d 4s to 3d 4s . 

(3) In addition to these, the bands tructure was also 
sensitive to whether the calculations were performed 
self -consistently or not. 

'The following potentials gave reasonably good band- 
structure results when compared with the previous theoretical 
and experimental results; 

71 ; constructed nonself-cons is tently with a 
com iguration 3d 4s and exchange parameter a = 2/3 [1]. 

72 ; constructed nonself -consistently with a configu- 
ration 3d'4s^ and exchange parameter a = 1. 
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Y3 : constructed self --consistently with a configura- 
tion 3d^4s and exchange parameter c: = 2/5 [2]. 

The Fermi surface constructea for these potentials 
were found to be in agreement with those constructed by 
Kattheiss [10], 'Jalroh [7] ant: pa paeons tantopoulos [2]. The 
comparison of the calculated ratios of the dimensions of 
the hole ellipsoid centred at the symmetry point IT agreed 
fairly well with the experimental values reported [11], 

Some of the previous calculations had reported the 
possibility of necks along the <110 > direction in the Fermi 
surface hole sheets connecting the JG- arms along <100 > direc- 
tions and the ellipsoids at IT [10], however, no such necks 
were obtained by us for any of the above three potentials. 

Calculations of the EKD were performed for the above 
3 potentials, detailed calculations being completed for 1/1 
and Y2. For these the EMU was calculated along 12 p direc- 
tions in the (100) and (110) planes. The EKE thus obtained 
was found to obey the symmetry selection rules proposed by 

_ _ 4 - 

Mijnarends [13]. The total EMD along the various p direc- 
tions could be accounted for in terms of the symmetry 
selection rules f Fermi surface effects and the intrinsic 
nature of the individual band wavef unctions . The behaviour 
of the EMI) was brought out through isodenses in the (100) 
ana (110) planes. It is suggested that the directional 
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Compton profiles from V single crystals could be experimentally- 
measured and unfolded to obtain the EMI) and then compared with 
the theoretical values, because of the sensitivity of the EMD 
to various features of the I : S (in particular the neck along 
the <110 > direction in the 3rd band hole sheet), this is 
especially desirable in vana.dium. 

Even though the bands true ture calculations are sensitive 
to the variations in the 3d population, the exchange parameter 
a, and whether the calculations are self-consistent or not, some 
of these variations have a cancellation with each other. Ihus 
for the potentials VI and V2, there are no significant differ- 
ences, while for VI and V3, some differences are observed along 
the < 110 > direction. 

It might appear that the results of our calculations of the 
EMD in vanadium are not so conclusive in distinguishing between 
the various potentials VI, V2 and V3. It should be pointed out, 
once again, that in going from VI to V2, the cancellation effect 
has taken place because the decrease in the d bandwidth due to 
the change in cc from 2/3 to 1 has almost compensated the increase 
in d bandwidth due to the 3d population going from 3 to 4. 
Similarity between the results of VI and V2 is, therefore, not 
surprising. Unfortunately, our attempts to calculate the EMD 
for some non-self consistent potentials with the exchange 
parameter a=2/3, e.g., 3d^4s^ a=2/3 could not succeed because 
of the computational limitations of Hubbard's approximate scheme. 
This limitation can be removed by using a more accurate method 
like the APW method to calculate the bandstructure and EMD. 
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Such an Ap¥ calculation using different crystal potentials can 
lead to a more conclusive choice amongst the potentials when 
compared with experimental results. This comparison for a metal 
like vanadium with a significant amis our opy in the EMD would be 
more meaningful if it involves a number of p directions. 

In our opinion,, the important conclusion brought out by the 
present work is as follows; For a metal like V, it is the combi- 
nation of the effects of the electron configuration, exchange 
parameter and self-consistency (or non-self -consistency) for the 
different crystal potentials used that determines the total 
nature of the bands tructure , the FS and the EMI/. It is therefore 
not enough to consider the effects of the electron configuration 
alone as was done by Paakkari et al. [6], These authors analysed 
the experimental data for the Compton profiles from some 3d 
transition metals in the framework of the UFA model and con- 
cluded that a (3d) n ~ 2 ( 4s) 2 type of electron configuration is 
more appropriate for metals like !i o-iid V. We feel that such 
an analysis is somewhat incomplete and the conclusions drawn 
are somewhat hasty. Our results indicate that such conclusions 
can be drawn only after accurate bands tructure calculations 
employing different sets of electron configuration, exchange 
parameter and self-consistency are made and results carefully 
compared with experimental data. 
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Chapter 5 


DIRECTIONAL COMPTON PROPILES ALL 
THEIR ANISOTROPIES FOR VANADIUM 


5 . 1 Introdu ction 

In the preceding Chapter we presented our theoretical 

results for the EELu in V in the (ICO) and (110) planes and 

discussed these results in terms of the bandstructure and 

Fermi surface of Y. It was pointed out that rhe EML curves 

obtained by us for the various p-directions could not be 

directly compared with the experimental results because no 

information about the directional EMD, p(p), unfolded from 

the experimental Compton profiles, (p ) , has beer- 

reported in the literature. There exists another meeting 

ground between theory and experiment and this is through a 

comparison of the directional profiles ( p ) and their 

anisotropies (or difference curves) Aj(p ) defined by 

z 


AJ(p 2 ) - J i 1 <] £ fiA(P z ) “ Jhhl^s^ 


(5.1) 


It has been demonstrated, particularly in the case of transi 
tion metals [1], that as compared to the single profiles the 



165 


anisotropy curves provide perhaps a more reliable and. 
sensitive method to compare the theory with experiment. The 
anisotropy curves are to te prefer3."'ea because of the fol- 
lowing two important a o. vantages ; (i) In the transition metals 

the core electron contribution ro the profile J(p ) is quite 

2 

large (about 58% at p =0 for V) . In a. single directional 
profile of a transition metal the. band profile is, therefore, 
masked by the core profile. However, the core profiles 
cancel out (in the approximation of the isotropic core EMi)) 
in the Aj(p ) curves and thus the anisotropy curves are more 
sensitive to the effects of band electrons. (ii) Secondly, 
the directional Compton profiles are affected by the multiple 
Compton scattering in the sample, the effect showing up as 
the reduction of the CP height, AJ(0), as well as the broaden- 
ing of the profile J(p ) itself. Although corrections for 
multiple scattering effects can be made these corrections 
are alwaj^s somewhat uncertain. If the sample thickness is 
held constant, the anisotropy curves AJ(p ) should not be 
significantly affected by the multiple scattering correc- 
tions. In view of the above considerations ve shall present 
in this Chapter the results of directional Compton profiles 

and their anisotropies Aj(p ) obtained from the Slid, p(p), 

2 

ca 1c ula te d by us . 



5 . 2 M ethod o f Calcu latio n 

The calculation of the directional Compton profile 
J-tviCP ) fz*om the theoretical Si'ir, p(p), requires an 

XjHi — L. 2 

evaluation of the double (or surface) integral 

+ » 

J hkP p z/ = ^ / p(p) dp x dp y (5.2) 

The integration involved in (13.2) requires a knowledge 
of p(p) at several points throughout the momentum space. 

This, in turn, involves a tremendous computational effort. 

An alternate but approximate method, suggested by Seth and 
Ellis [2] consists of expanding p(p) as a series of lattice 
harmonics. The outline of this method is given below. 

It was first suggested by Fdj'narenas [1,3] that the 
EMi p(p) can be expanded as 

P(P) = l '^..(p) -iv ^ Q tr (5.3) 

l,v xv p 

where F (fi->) are the lattice harmonics which form an 
1 v P 

orthonormal set of invariant linear combinations of spherical 

harmonics of order 1 while denotes the orientation of 

P 

p. .different harmonics of the same order ana symmetry type 
are distinguished by the label v . In the case of cubic 
lattices (e.g. bcc Y) only the harmonics having 1 = 0,4,6,8,10, 
12 (twice), 14,... will contribute. 

Next two orthogonal coordinate systems A and R are 

C ct 

defined. A system fixed to the crystal ana denoted by 
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H = (£ 5 t ) 5 S ) is used as the reference system for defining 
the lattice harmonics.' Another coordinate system R & =(x,y,z) 
is assumed to be fixed to the apparatus so that the z-axis 
is along the scattering vector, i. rotation of the coordinate 
system R through the .drier angles a, j3 and y can make the 
two systems £ w and Ii c coincide. fie angle a is a measure of 
the rotation of the sample around the scattering vector 
while the angles ( p ; y) denote the polar and azimuthal 
angles of the z-axis (i.e. the scattering vector). 


rhe lattice harmonics can be further expanded 

in terms of the spherical harmonics (ft-*) through the 

lm p ° 

relation [3] 


\v la f> = l ^ 


(5.4) 


It has been shown by Lijnarenas [3] that the directional 
Compton profile Jp a (p z ) can be written as 


V (p 4 = z% l 


l,v 


? lv (|3 5 a) S lv (p z ) 


where 


co 


~ , ^ , p iv^ p ^ k-jCpg/p) p dp 


IP 


r l^z' 


z 


(5.5) 


(5.6) 


where P^'s are the Legendre polynomials. 

In his original method Mi jna rands [3j used the scheme 
of Eq. (5.3)”(5.6) to unfold p(p) from a set of experimentally 
measured Jp.„(p z ). he followed the scheme of Seth and Ellis [2] 


9 
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and Panhere [4] to obtain J o„(p„) from p(p). This method 
L J pcc z 

cons is tea of the following steps . 

i) The starting point of our calculations was the set of 
p(p) distributions along the 12 p-directions (Pig. 
4-24) obtained by us for the VI type of the crystal 
potential of V. Ph.es e p(p) curves are shown in 
Pigs. 4.11, 4.16, 4.21 and 4.25-4.33 (Chapter 4). 

ii) As a next step we used hq. (5.3) to expand p(p) in 

terms of p, (p) and P_ (ft-*.) to obtain a set of linear 
lv ~ ±v p 

equations. The expansion (5.3) was truncated at 

1 = 12 . 


iii) This set of linear equations was next solved by using 
a least squares fit and an inversion scheme to obtain 



iv) 'The set of p, (p) functions so obtained was next used 

lv 

in So, (5.6) to construct V (P Z ) functions which 
through 5q. (5.5) lea to the construction of the 


direct io nal prof il e s 


J- (p ) 

•- (y v u rr ' 

(L (X 


5 . 3 Results and_ b is c uss ion 

The g and P functions obtained by the above 

1 V JLV 

calculations are shown in Pigs. 5.1, 5.2 and 5.3 for 1=0, 

4 and 6. The function g^ describes the spherically averaged 
part of the Compton profile (see Eq. 5.6) while the higher 




Cass VI 



F±S * 5,2 * Th ® Kubi0 harmoniG components g 4 (p) and gjp) of the Compton 
profiler {case ¥1), 



Case VI 
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(1 >_ 4) g-, -functions indicate the non- 
anise tropic part of the Compton pro ail 
function describes the spherically ave: 
while the higher (1 >_ 4) p j functions 


spherical or the 
e. Similarly, the p^ 
x^ed part of the EME 
describe the anisotro- 


pic part of the EEtu. 


Examination of i'i^. 5-2 reveals that ana gg are 
of comparable magnituue and they augment each other in the 
region p < 1.2 a.u. but in the higher momentum region they 
cancel out. the curve for Pg shows more discontinuous 
behaviour than that for P„ . i'ne structure shown by P r 
arises from the spikes in the Elio along, the < 110 > direction 
and is related to the Fermi surface topology (e.g. dipping 
oi the 2^ band below the Fermi level) ana its Umklapp images. 

She results for g^ ana p 1 functions shown in Figs. 5.1 
5.3 have one usefulness. ihe unfolding of the experimental 
Compton profiles by the lattice harmonic expansion method 
[1,5] also yields the plots of g^ and curves. Examina- 
tion of plots similar to those shown in Figs. 5. 1-5. 5 can, 

therefore, help design the Compton profile experiment by 

. . 

choosing the pi-directions (Bq. 5.5) for CP measurements. 

ihe Compton profiles from polycrystal- 
line samples of Y has been measure a by hanninen and paakkari 
[5] and were later corrected for the multiple scattering 
effects [6 j . In order to compare our theoretical profile 
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■Table 5.1* 

Comparison of the 
profile (present 
mental values for 

spherically averaged Compton 
calculations) and the experi- 
■ poly crystalline V of he f. [6] 

P z in a.u. 


^sph^ -^z^ 
present 

J s ph ( Pz } 

calc. Ref. [6] 

0 . 00 


5.629 

5.405. +.0.08 

0. 05 


5.630 

- 

0.10 


5.607 

5.387 

0.13 


5.575 

- 

0.20 


5.529 

5.322 

0.30 


5.403 

5.210 

0.40 


'5.224 

5.055 

0.50 


4.999 

4.861 

0.50 


4.728 

4 . 636 

0.70 


4.434 

4 . 386 

0.80 


4.135 • 

4.121 

0.90 


- 

3.849 

1.00 


3.518 

3.576 + 0.07 

1.10 


- 

3.308 

1.20 


2c 904 

3.049 

1.30 


- 

2.802 

1.40 


2 o 48^J 

2.572 

H 

* 

VJ) 

c 


“ 

2.360 

1.60 


2o 023 

2.170 

1.70 


- 

2.002 

1 „ GO 


1.710 

1.855 

1.90 


- 

1.730 

2.00 


1.482 

1.623 + 0.05 

2.20 


- 

1.452 

2.40 


1.128 

1.309 
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with, their results we proceeded as follows. Using the 
spherically averaged band profile [given by gQ(p z ) of 
Pig. 5.1] we first normalised it. so that 

2.4 a.u. 

j S 0 (P z ) d P z = 2.2774 electrons 

o 

where the number 2.2774 electrons was obtained by finding 

*5 2 

the area under the free-atom (3d) (4s) profile of Biggs 

et al. [7] between 0 and 2.4 a.u. To this normalised band 

profile we adaea the 'core' profile due to the (is) (2s) 

(2p) ^(3s) ^(3’p)^ electrons from the iables of Biggs et al.[7j. 

Such a total profile calculated for polycrystalline V is 

shown In Big. 5.4 which also shows the experimental data 

by paakkari et al. [6]. Both these profiles are normalised 

to the same area between 0 and 2.4 a.u. and are also compared 

in liable 5.1. The theoretical profile was convoluted with 

the residual instrument function of the experiment [6], It 

is observed from Pig. 5.4 that our theoretical calculations 

slightly overestimate the contribution at low (p < 0.6 a.u.) 

z 

momenta but underestimates it at the high (p 2 > 1*1 a.u.) 
momenta. It may be pointed out that the theoretical Compton 
profile of V calculated by using the renormalised -free-atom 
model [6] also gave lower values in the momentum region 
P z = 0.5-1. 5 a.u. The slight disagreement between our 
calculated results and the experiment could be due to the 
inadequacy of the independent particle model used by us or 
due to some other reasons which are suggested below. 
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In the calculation of the spherically averaged as well 
as the directional profiles by the lattice harmonic expansion 
method we had made two approximations, firstly the p(p) 
distributions used by us were calculated for the 12 ^-direc- 
tions (life. 4.24) all of which were lying in the symmetry 
planes, the (100) and the (110) planes. Thus the p(p) 
distributions along the general p-directions (those removed 
from the (100) and (110) planes) were not considered while 
constructing the spherically averaged profile ana this may 
lead to an incomplete spherical average. Secondly the 
accuracy in the calculation of the EMu by the Hubbard - 
rli jnarends method is limited at high momenta (p > 3 a.u. ) 
and we had therefore calculated the p(p) values upto 
p v3 a.u. As a result the integrations involved in Eq. (5.6) 
were necessarily truncated at T? max ^3 a .n. and this may 
lean to some error. In addition, another source of error 
could be the muffin-tin approximation to the crystal poten- 
tial. It has been pointed out [0,S] that the non-constancy 
of the potential in the interstitial region ('warped muffin- 
tin’) and the non-sphericity of the potential inside the 
muffin-tin sphere can affect the band structure results. 
Although no estimate of these effects on the EMu are 
available, we feel that some of the disagreement observed 
in Fig. 5.4 could also arise partly from these effects. 
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vv'akoh and Yamashita [10] have calculated the spherically 
averaged Compton profile of Y hy a direct integration involved 
in Eq. (5.2) and obtained a good agreement with the earlier 
results of Paakkari et al. [see Ref, 10] which were perhaps 
not corrected for multiple scattering. It would be interest- 
ing to know whether a comparison of their theoretical profile 
xdth the re, vised experimental data [6] corrected for multiple 
scattering shows a behaviour similar to Pig, 5.4« The 
directional profiles obtained by us for the principal symmetry 
directions <100>,<110> and <111 > are shown in Rig. 5.5 and 
their general features agree with the theoretical profiles 
calculated by vlakoh et al. [11] employing the API'.' method and 
by Callaway et al. [12] employing the LOGO method. 


The most striking feature of the directional Compton 

profiles of Y is the dip shown by the J]_qq(p 2 ) profile at 

p < 0.3 a.u. (Pig. 5.5). i‘o understand the origin of this 

dip it is not enough to examine the p(p) curves (Chapter 4). 

The directional profiles a,re constructed by the integration 

(Eq. 5.2) in the momentum space and one should therefore 

examine the EMI, p(p>), in different (p„ , p ) planes of Integra-- 

tion (normal to the hkl direction) for p = constant. However. 

z 

our present results of the EICu are restricted to the (100) 
and (110) planes, following Vakoh. and Yamashita [10], a fair 
understanding can be obtained by examining the Fermi surface 




The ache mat ic diagram of the Pd topology 
showing the various symmetry planes [10], 


Pig. 5.t> : 
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topology of V in the repeated zone scheme (Fig* 5.6)* It 
must be remembered, however, that in addition to the Fermi 
surface topology, the effect of the symmetry selection rule 
(explained in § 4.4) "trill have to be taken into account for a 
complete understanding of the problem* 


The concave structure (or the dip) shown by the 
profile at p z <\,0.0 a.u. and again at p^^l.l a.u. arises pri- 
marily because of the jungle-gym structure of the Fermi surface 
of V (Chapter 2). In part 2 of the Fig. 5.5 are shown the 
planes of integration for the scattering vector along the 
[100] direction. Plane a passes through the point r(000) or 

p = 0 while the plane a' passes through the point H (100) or 
z - 

2 IT 

P z = g— = 1.034 a.u. '.These planes a and a ! contain the jungle- 
gym hole arms (or tubes) in two directions as well as the 
ellirosoidal holes around N. These unoccupied states will not 
contribute to the SMD in these planes. As one moves the plane 
of integration from a towards a', say at p v 0.3 a.u. the 

Z 


integration planes do not contain the jungle-gym hole arms 

(or the hole ellipsoids) any more. There is, therefore, a 

greater contribution to the EMD in this plane (as compared to 

a) and as a consequence the Compton profile at jp^^O.5 a.u. 

shows a value which is higher than at p = 0.0 a.u. This 

explains the concave part in the J, rn (p ) below p = 0.3 a.u. 

XUu z z 

The Umklapp image of this concave part occurs around p =1.1 a.u. 

z 

as another concave curvature . 
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Similar considerations for the J-^qCp ) show that the 
two planes b and b ' (Pig. 5.6) contain the jungle-gym hole 


arms in one 
profile 
concave part 
= 0,774 a.u. 


(ana not two) direction only. As a result the 

(p ) displays a less pronounced (or shallow) 
z 

at p =0 ana again weakly at p = ( 2Tc/a) /^"2 

3 3 s ’ 

through the Umklapp process. Compared, to the 


ana P r °£il es the profile shows no dips or 

concave structure because the planes c ana o' (Fig. 5-6) 


normal to the [111] direction do not contain the jungle-gym 


hole arms. 


When compared to the directional profiles (Fig. 5.7) 
calculated by Uakoh et al. [11] the present results (Fig. 5-5) 
show a qualitative agreement. However, the theoretical 
profiles of Uakoh et ai. [ll] give somewhat lower values at 
low momenta and slightly higher values at high momenta. 
Unfortunately, there are no Fables of experimentally measured 
directional profiles of Y reported in the literature and 
hence a direct comparison between experimental and theoretical 
profiles obtained by Wakoh et al. [11] or by us is not 
possible. VJe feel that the differences between the results 
of Wakoh et al. [11] and the present work could be due to 
the different methods followed by wakoh et al. [11] to 
construct the crystal potential (they have constructed the 
potential phenomenologically fitting state -dependent poten- 
tials to experimental FS data) , to calculate the bandstructure 





Tho d Irociional Uov too. profiles calculated b 
!A: ' on a% •'• ,1 , j. II j . 






J 100 

Present results (VI ) 



or- 









Tne 
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and EELu (they have used the ApV. 7 method) ana to construct 
J hkl( p z) prom p (p)* ' Bie l ast feature could he most important 
because Wakoh et al. [11] have used a- direct integration in 
the momentum space (2a. 5.2) while we have used the approxi- 
mate method of lattice harmonic expansion. 

As pointed out earlier the experimental data for the 
directional profiles for V was reporteu by ferasaki et al. [15] 
in the form of anisotropy or difference curves and is shown 
in Fig. 5.8a. The corresponding anisotropy curves obtained 
from our directiona.1 profiles of Fig. 5-5 are shown in Figs. 
5.9? 5.10 ana 5.11. 


A quantitative comparison between the experimental and 
the present theoretical anisotropy curves is not possible as 
Terasaki et al. [15] have not reportea any numerical tables 
in their paper. However, a qualitative comparison of the 
results in Fig. 5.8a and Figs. 5.9-5.11 shows a satisfactory 
agreement except for the ( J-^n" ^100^ * curv ’ es i n the high 
(p 2 >1. a.u. ) momentum region. The amplitudes of the sub- 
maxima in the experimental and theoretical anisotropy curves 
show a difference which could be ascribed to the effect of 


smoothing by the detector resolution function. The anisotropy 
curves calculated by us show very good agreement with the 


calculations of Wakoh et al. 
curves at very low momenta. 
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curves (Ilg, 5.8a) show a better agreement with our results. 
It may be pointed out that neither the experimental curves 
(Pig. 5.8a) nor the present theoretical curves (Pig. 5.9- 
5.11) agree with atomic model of Weiss [14] whose results 
(Pig. 5.8b) did not predict any submaxima and subminima in 
the anisotropy curves (Pig. 5.7b). 
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COMPTON PROS' ILLS OP MALI Till LBUESRILS 


f . 1 Introduction 

Metal hydrides (deuterides) are solid compounds of 
hydrogen (deuterium) with metals. In recent years the study 
of metal hydrides has evoked great interest because of their 
possible use as storage elements of hydrogen for both fixed 
plant and vehicular applications [1], In fusion reactors 
vanadium is considered as a construction material and severe 
constructional problems can arise due to the solution of 
hydrogen in vanadium [2]. In view of these technological 
applications, the physical and other properties of metal 
hydrides are being increasingly studied by different experi- 
mental techniques. 

A fundamental and important question in these studies 
is the electronic structure of the metal hydrides. Although 
other experimental methods like photoemission studies, soft 
X-rasr emission spectroscopy etc. are useful to investigate 
the electronic structure of metal hydrides, it has recently 
been demonstrated [ 3 - 1 ] that the Compton scattering method 
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provides a sensitive prole to study this aspect. She advan- 
tage of Compton profiles is that they are more sensitive to 
the electronic states ox the conduction electrons but much 
less sensitive to the lattice distortions introduced during 
the process of loading hydrogen into the metals. 

2w o models have been proposed in the literature [8] 
for the electronic state assumed by hydrogen in the metal 
hydride. In the protonic model the hydrogen (or deuterium) 
atom is assumed to donate its electron to the conduction 
band of the host metal while in the anionic model the 
hydrogen is assumed to remove an electron from the conduction 
band to form an anion in a bound Is state. Ab initio band- 
structure calculations [9-12] of palladium hydrides show that 
the electronic structure of metal hydrides is more complicated 
than either of the above two simple models. 

Previous Compton profile measurements on polycrystal- 
line YEq ^ [3] and on single crystalline Nbijg g [4], 

IHdHq yg [5] and ITbHgj [6] supported the protonic model for 
these hydrides. Recent studies [7] of the Compton profiles 
of polycrystalline P&Hq ^ and VHq ^ support a 

model which takes into account the creation of metal hydrogen 
bonding states below the Fermi level by the introduction of 
hydrogen in the host lattice. 

Among the transition metal hydrides, the vanadium 
(or niobium) hydrides are expected to show interesting 
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bancistructure changes because of the complicated bandstructure 
and Fermi surface of the host metal. flic discussion in the 
previous Chapters has shown that the- electron momentum 
distribution (ML) in vanadium is very sensitive to the 
electron states near the Fermi level. Tic, therefore , thought 
it interesting to extend our calculations of the FML for pure 
vanadium to vanadium deuteride, VLq ±To ab in itio band- 

structure calculations of the EILD in vanadium hyaride (or 
deuteride) are reported in literature yet. Gur present 
calculations can provide a theoretical basis to compare the 
experimental Compton profiles from the single crystalline 
vanaciium deuteride (hydride). The present calculations are 
particularly applicable to vanadium deuteride -VLq ^ 
because for this system (at the room temperature) the vanadium 
lattice does not undergo a phase transition by deuterium 
loading- the other hand, the phase diagram of VH is 

■vorj complicated [13j ana it indicates that at room tempera- 
ture pure V shows the symmetry BCC (cs-phaso) while VH-. 

VJ « / 

shows a possible symmetry 3CT (e-phase) [14]. 

6.2 Meth od of Calcula tion 

Perhaps the most correct approach to the calculation 
of the Compton profile of a metal hydride is to perform an 
appropriate bandstructure calculation of the hydride itself 
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and calculate the EMD from, the band wavef unctions so 
obtained. The Compton profiles could then be obtained from 
the double integration in the k -space. However the band- 
structure calculations of the metal hyarides performed so 
far [9-12] ao not report the wavef unctions or the Ml. 

The basic iaea behind this calculation is taken from 
the suggestion made by Kanhere ana bingru [15]. These 
authors stuaieu the systematics of the theoretical electron 
momentum aistribut ions in the 3d transition metals V, Cr, 

Pe, hi ana Cu. They pointed out that along a. particular 
< hkl > direction the nature of the EMD arising out of indi- 
vidual bands or the p. (p) curves (j being the band index) 

J 

does not change much in going from one metal to the next 
higher metal (i.e. from 7 to Or) having same (i.e. BCC or 
PCC) crystal structure. The total EMD, 

0CG -> 

p(p) = I p.(p) (6.1) 

3=1 5 

however, shows significant changes which can be correlated 
with the changes in S^. This interesting and systematic 
behaviour of the EMI) is similar to the behaviour of the 
density of states, n(E) , in the 3d metals. The systematics 
of the n(E) curves have given rise to the so-called ’rigid 
band model’ [16] which assumes a universal density of 
states, n(E) , for a given family of alloys. The Eermi 
energy (Bp) of the alloy is then defined through the relation 
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I'T = J ~ n(]j) dS (6.2) 

“ CO 

where IT is the number of conduction electrons in the alloy. 
:;e have made use of our results of the Eku for pure V to - 
calculate the EI-il ana Compton profile of Vuq ^ nn the basis 
of such a rigid band model. It is realised that an ultimate 
answer to the solution of this problem will have to come 
from the bandstructure calculation of the metal hydride 
system rather than an application of a rigid band approxima- 
tion. 

Earlier bandstructure calculations, particularly for 
palladium hydrides [9-12] have shown that the d-states in the 
host metal are almost unaffected by the hydrogen but the 
states possessing s-lihe character about the hydrogen will 
be modified. ‘Thus hydrogen- (or deuterium-) vanadium bonding 
states are created below the Permi level by the introduction 
of hydrogen (or deuterium) in the host lattice. 

To examine this effect we have calculated the Compton 
profiles and anisotropies for ^ by two different 

methods with the protonic model anc the rigid band approxima- 
tion as the basic features of the calculations. 

(a) In the first method we have assumed that in V1 q ^ 
deuterium atom donates all its electrons to the conduction 
band of vanadium. This we shall call as the simple protonic 
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model (SPK). fhe ener^r bands being tah en 
Fermi energy is raised from its value If. ^ or P 111 "® 
tial 71) to a higher value E-^, such that 


5.0 = J 11 (E) dE 


as rigid, t3ae 


Y (P° 


ten" 


(6.5a) 


"F 

5 .72 - / n(E) dE 


(6.5^) 


The total EMB p(p) was then calculeted sim:n: u 

^ (for 

Bq. (5.1) up to the new E^ using the same Pj(f) V °' 1U 

inciividual bands) that were obtained for pure v--.-c.di 

1 n "h P (5. 'f'o 

(potential 71). 'file total EMD was thus calc urn a. 

12 indirections in the (100) ana (110) planes shown a 


the 


Big. 4.24< 


’O" 


of electrons 


l • f ie d P r 1 

(b) In the second method to be called as the mo 1 
tonic model (MHi) we assumed that only a f rac cion of 

• ,-o the conduction 

electrons of the deuterium atom are donatea bu 

hands of vanadium while the remaining fraction 

are accommodated in newly created center ium-vanaaiuia 

, . rql h aS 

states formed below the Permi level. Switenaic-^ l-'J 
performed a bancs tructuro calculation on vana 
which supports the formation of such states anu Sii ° ws 
such vanadium-hydrogen bonding (or hybrid.) staoew Cc.n 


•ponding 


dium hydti^ 6 
that 


accommocate about 0.25 electrons per H atom. Using 


thi s 
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data for our MH'i we have assumed that 0.25 x 0.72 or 0.18 
electrons from the deuterium atom in Y'u r no go into the 
formation of vanadium-deuterium st?tes having Is character 
about the deuterium site . Their wavef unctions and EMf> are 
aescribed by free hyarogen atom wave! unctions . The remain- 
ing 0.75 x 0.72 or 0.54 electrons from the deuterium atom 
in TDq rj 2 are accommodated in the conauction band of 
vanadium raising its Fermi energy to a new value Ep such 
that 

Tp! * 

^F 

5.54 = / n(E) dE 

— co 

The EMI for these 5-54 electrons is calculated by s ummi ng 
Eq. (6.1) up to Ep. 

In both the cases (a) and (b) the directional Compton 
profiles (p ) were calculated from the total EMI), p(p), 

by following the method of Subic harmonic expansion for the 
12p-directions as described in Chapter 5. In the case (b) 
or MPM, the contribution to the Compton profile from the 
0.18 electrons occupying the vanadium- deuterium bonded 
states was obtained by using the tables of the Compton 
profiles for free hydrogen atom by Biggs et al. (17). 

The normalization of the area under the Compton profiles 
was done as follows. Firstly, the Compton profile of the 
free atom (3d)^(4s) valence electrons was taken from the 
calculations of Biggs et al. [17]. The area under this 



profile in the momentum region p = 0.0 to 2.4 a.u. was 
founa to be 
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2.4 a.Uc 


J 


iree 


atom 


(p z > 


d] 


IT = 2.2774 electrons. 


Next the Compton profiles * i.e., £ 0 (p z ) > ^hVl^z^ °alculated. 
by us for pure V metal (ce.se VI) were normalized to the 
same quantity If = 2.2774 electrons as follows, 


2.4 a. u. 

J J pure V d P 2 ~ N o 


For VDq rj 2 the Compton profiles obtainea by the simple 
protonic model (SPM) were normalized to a quantity 
N = 5.72 (N /5.0), so that 

Si o 


2.4 

/ 


o 


a.u. 


4? (p 2 ) 

“^ 0 .72 


a P z = »a 


For the modified protonic model (MEM) the normalizations of 
the Compton profiles for VFq corresponding to the conduc- 
tion band electrons (which in this case are 5.54 as explained 
before) and the 0.18 electrons from the deuterium atom 
(which form the s-like bonding states) are done separately, 
as explained below. 


For the 5.54 conduction or band electrons we put 


/' 4 8 " U ' « = 5-54 ( 5 o / 5 . 0 ) = N b . 



198 


'Do normalise tiie Compton profile clue to the C.18 electrons 
we started with the Compton profile for 'die hydrogen atom, 
as given by Biggs et : ~1. [17] end obtained a number 

IT i r om 

2.4 a.u. 

f Jt-(p ) do = :T = O.dG .33 electrons. 

j hi. x z - 2 C 


'Hie normalization for the C.18 electrons of the deuterium 
atom which form the s-like bonding states in YD P r , 0 in the 
moaified protonic model was then carried out as follows 


2.4 a.u. 

/ 


T MBi /_ x d _ 
J 0.18D p z dp z 


Finally the total profile for VI 
model was obtained as 


= 0.18 x IT 

C- 

q j 2 in bhe modified protonic 


^0.72 


(P z ) 


J- 


MPM 


Vb, 


0.54 


(pJ 


+ 


J-IHi 
d 0. 181) 


W z > 


It may be pointed out that in the present calculations there 
is basically no difference between hydrides (YH ) and deute- 

-A. 

rides (YD ) . However, our theoretical results are more 
directly applicable to the deuteride YDq ^2 rather than to 
VHq Y 2 because of the fact that the former system does not 
undergo a phase transition from pure Y while the latter 
system does . 

It is also important to point out the difference between 
the method of present calculation and that of Lasser and 
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Lengeler [7]. Although these authors describe their calcula- 
tions as a bands tructure calculations, the basis of their 
work is the Compton profile theoretically calculated by 
¥akoh et al . [18] for pure V using the API/ banastructure 
method. L'dsser and Lengeler thus consider the Compton pro- 
files as ’rigid' and adjust the Compton profile contributions 
according to the different models used by them. On the 
other hand, we feel that in a metal like V where the EMI) is 
so sensitive to the changes in the banastructure and Fermi 
energy one should treat EMI) as more basic than the Compton 
profiles. Accordingly we have first calculated the EMI in 
YL 0 ^2 as modified by the transfer of electrons from 
deuterium to the host metal, ihe changes in the Compton 
profiles in going from V to VLq ^ are then calculated by 
us from the modified EKD under the two models (i.e. Sffi and 
MEM). Further our calculations, which start from the total 
EMI) along various ^"-directions, lead to the prediction of 
directional Compton profiles as well as anisotropies for the 
single crystals. The present calculations thus provide a 
wider basis for a comparison between theory and experiment 
and allow a more sensitive test of the two models (i.e. 

SPM and MEM) proposed. 

6.3 Results and Discussion 

The results of our calculations for the polycrystalline 
YLq rj 2 are shown in Fig. 6.1 in the form of a difference curve 
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aJ(p 2 ) = J VD (p a ) - J v (p 2 )„ The isotropic profiles 
(corresponding to the poly crystalline Vn,-, or V) were 
obtained from the respective g Q functions (Do. 5.6) and were 
normalised as explained in the previous section, fhe theore- 
tical profiles have not been convoluted by the energy resolu- 
tion function of the detector system. i'he curve narked SIM 
refers to the simple protonic model calculation while the 
other curve marked MPIi describes the modified protonic model 
calculation. There is no experimental measurement of the 
Compton profile for VI>q reported in the literature, 
nevertheless we have obtained the Compton profile for 
VIq ^2 by using the experimental data of Lasser and lengeler 
for VI)q rjrj [7] and normalising the area appropriately. The 
experimental data points for VDq so obtained are shown 
by (.) marks in Fig. 6.1. Earlier work has already shown 
that neither the atomic model [7] nor the anionic model 
[5-7] is able to describe the Compton profiles of transition 
metal hydrides satisfactorily. A cor peris on between our 
theoretical curves and experimental values (fig. 6.1) confirm 
the earlier finding [7] that a simple protonic model cannot 
correctly describe the difference curves AJ(p_) in vanadium. 
Our present results along with those of L&sser and Lengeler 
[7] suggest that a modified protonic model which takes into 
account the formation of hydrogen-vanadium bonding states 
below the Fermi level is more successful in describing the 





100 



‘i-*T f 1 O r > Q T .i O Jf 
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observed experimental results. Hie residual disagreement 
between the experiment ana MPEa values (Pig. 6.1) could be 
due to the approximations made by us in the PPM and those 
involved in deriving J^p^) from the 311 (Chapter 5) . 

It was seen in Chapter 5 that in vanadium, the directio- 
nal Compton profile along the <100 > direction reflected the 
complicated Permi surface topology of vanadium - particularly 
the features like the jungLa-gym hole tubes along the <100 > 
directions, lie have therefore compared the theoretical 
directional profiles J^OO^z^ f° r ^ a ad Vjjq ^ (SPEI) in 
Pig. 6.2. It is seen that the concavity of J-)qq(p z ) in 

the p range from 0.0 to 0.5 a.u. observed in vanadium 
z 

becomes more flat in Yl) r n 0 . ‘Phis effect is also seen in 

u . / a Y'j) 

the difference profiles AJ 1 qq(p z ) = ^00°^^ ” ^100 

Pig. 6.5. Phis behaviour is easily understood in terms of 

the bands true ture and Permi surface of vanadium. 


In the construction of the Compton profile (CP) J^OO^z^ 

the integration planes of the PHD (3q. 5-2) are parallel to 

the (100) plane. Por p = 0, the (ICO) integration plane 

z 

passing through r contains* in the case of vanadium, four 
jungle-gym arms (due to the third band Permi surface) along 
the and +p^. directions. Phe density of occupied states 
for this plane is relatively lower and the comes ponding 
EMI has the features shown in Pig. 4.54. As a result, the 


9 



205 


Jl^o(P 7 ) curve shows a concave part in the region p z =0.0-0.3 
a.u. In V'Vq ^ 2 ' however, the electrons from the deuterium 
atom fill up the third bana considerably thus leading to a 
shrinkage of jungle -gym hole arms. Thus, in the (100) plane 
of Vj jq rj 2 passing through p z = 0, there are relatively more 
occupied states and the contribution to the EKD is corres- 
ponding!;?' larger. This results in a more flat behaviour of 
the J-|_qq curve for VDq ^ region *p z ; < 0.3 a.u. 

Similarly at the points H' and the higher r.-point along the 
<100 > direction, we have a more flat behaviour in the Compton 
profile for VDq y ^ compared to that in vanadium. Further, 
the comparison of the JpQ()(P z ) profiles in Nb and FbHQ yg [5] 
show similar features which are explained by arguments same 
as above. This is not surprising in view of the close 
similarity between the Fermi surfaces (and b anas true ture) 
of Y and Nb. 


The difference curves 


vjj n 79 r 

A Jll 0 (P z ) = J no* ( p z ^ “ J 110^ p z^ 


AJ lll^ p z^ ” J lll' 72 ^z^ J lll^ P z' > 


' 111 *' *z' 


are shown in Figs. 6.4 and 6.5. These results show that 
compared to the difference curve AJ-^y, the curves AJ^OO 
and A Jy 2 _o s* 10 ™’ more significant changes in going from the 




( z d) 0Ll rv 









0.03 





0.03 



The theoretical Compton profile anisotropy curve 
(SEW and MPM) for ¥D Q 
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simple protonic model to the modified protonic model. The 

anisotropy curves ( ^hhl”^h ' k *1 1 ^ ^ or 72 '^ e "^ wo 

ferent models (Sffi and I4PM) are given in Pigs. 6. 6-6. 8. Here 
we find that the two models predict a significant difference 
in the anisotropy curve Thus it appears that the 

comparison of the theoretical results for the difference 
curves a J 1 qq» A j ^ as vel - 33 anisotropy curve ( Jppo“^lll^ 

with experimental values for these canes can throw more light 
on the relative merits of the two models. 


No measurements of the directional Compton profiles of 
Vl (or VH ) have been reported yet. Such measurements on 

X XL 4 ~ l ” 

NbH and Nbjj. [3-6] have already added to our understanding 
of the electronic structure of these hydrides. He suggest 
that the measurement oi the Compton profiles from the single 
crystals of VS and VH be undertaken. A comparison of these 
experimental results with the results of our calculations, 
i.e.. Pigs. 6. 3-6. 8 can then provide a basis to understand 
more clearly the electronic states in Vcaiadium neuter ides 
and hydrides . 


Summary 

Compton profile measurements have been found to be a 
suitable probe to study the electronic structure of metal 
hydrides and deuterides. Sven though earlier experimental 
studies favoured the protonic model, recent experiments [7] 
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have shown, that this monel cannot correctly explain the 
electronic band structure of metal hydrides. In the case of 
some transition metals, bands tructure calculations have shown 
that only a fraction of the hyarogen (or deuterium) atom 
electrons are transferred to the conduction band while the 
rest of the electrons form Is bonding states with the 
electrons below the Permi surface of the host metal. 

"vie have used our results for pure vanadium to investi- 
gate this effect. Compton profile calculations for VIq ^ 
were done by two different methods -under the rigid band 
approximation (i.e., keeping the partial momentum densities 
rigid) and with the protonic model as the basis. In one, 
all the electrons of the deuterium atom are transferred to 
the conduction band (SPM) while in the other, only a fraction 
are transferred, the rest being assumed to form Is bonding 
states with the electrons below the ?e rmi surface of the 
host metal (MPK) . She results of the calculation show that 
the I-lII'i is more successful, in explaining the experimental 
values obtained for pclycrysta'line vanadium deuteride. ihe 

difference curves i J = J Tn - J, r as well a.s the anisotropies 

“^ 0.72 V 

^"hkl “ ^h'k'l 1 ^ or ^ ne aEL ^ s k°“' <r significant differ- 
ences. It is suggested that comparison of these theoretical 
results with future experimental results on single --crystalline 
Tjj q rj2 could throw more light on the electronic structure of 
hydrides ana deuterides. 
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